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ABSTRACT 


Capillary  forces  may  have  a  major  effect  on  the  behaviour 
of  immiscible  displacement,  because  they  affect  the  saturation 
profiles  and,  consequently,  the  ultimate  oil  recovery.  A  quick 
method  for  capillary  pressure  measurement  was  investigated.  It  makes 
use  of  a  centrifuge,  but  differs  from  the  conventional  centrifuge 
method  in  that,  instead  of  using  small  core  samples,  measurements 
are  made  on  large  samples  cut  into  slices.  It  was,  therefore, 
possible  to  obtain  simultaneously  several  points  on  the  capillary 
pressure  -  average  saturation  curves.  These  data  do  not,  however, 
represent  the  true  capillary  properties  of  the  system  under  study. 
That  is,  they  are  affected  by  several  experimental  parameters,  such 
as  size  of  the  slices  and  whether  one  or  several  speeds  of  rotation 
are  used  for  obtaining  a  curve.  These  effects  were  studied  in  an 
attempt  to  get  a  better  understanding  of  what  the  experimental  data 
thus  obtained  really  represent.  Finally,  a  new  mathematical  model 
defined  by  four  parameters  is  proposed  to  fit  the  data.  It  was 
shown  that  a  reasonable  estimate  of  the  area  under  the  curves  may 
thus  be  obtained.  This  term  is  strongly  related  to  the  parameter,  a, 
used  in  the  definition  of  the  capillary  number  which  appears  in  the 
immiscible  displacement  and  fractional  flow  equations. 

Unsuccessful  attempts  were  made  to  obtain  imbibition 
capillary  pressure  data.  These  failures  seem  to  be  the  consequence 
of  limitations  in  the  equipment  used. 
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CHAPTER  I 


INTRODUCTION 


Two  or  three  immiscible  fluids  are  generally  present  in 
petroleum  reservoirs.  As  a  consequence,  forces  are  created  because 
of  differences  in  the  thermodynamic  properties,  such  as  chemical 
potentials,  associated  with  the  various  components  present  in  the 
phases.  These  differences  give  rise  to  a  physical  interface,  which 
has  a  finite  thickness,  and  within  which  the  thermodynamic  properties 
are  changing  continuously.  The  forces  created  at  the  interfaces 
define  the  capillary  properties  of  the  system. 

The  initial  distribution  of  the  fluids  present  in  a 
natural  reservoir  depends  on  the  capillary  properties  of  the  rock 
and  fluids  system.  Although  the  magnitude  of  the  capillary  forces 
is  generally  small,  these  forces  may  have,  in  some  instances,  a  major 
effect  on  the  behaviour  of  immiscible  displacement;  the  saturation 
profiles,  which  determine  the  ultimate  oil  recovery,  are  dependent 
upon  the  capillary  effects.  It  is,  therefore,  necessary  to  know  the 
capillary  pressure-saturation  relationship  for  every  system  under 
study.  This  involves  the  acquisition  of  accurate  experimental  data, 
as  well  as  a  good  understanding  of  what  these  data  really  represent. 

Leverett  (1)  used  the  definition  of  capillary  pressure, 
together  with  Darcy's  law  written  for  two  phases  and  the  continuity 
equation,  to  derive  his  fractional  flow  equation.  In  a  subsequent 
paper,  Buckley  and  Leverett  (2)  solved  this  equation,  assuming 
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that  the  effects  due  to  gravity  and  capillary  pressure  were 
negligible.  The  solution  of  this  problem  led  to  a  saturation  profile 
which  was  triple-valued  over  a  portion  of  its  length.  Buckley  and 
Leverett  (3)  recognized  the  physical  impossibility  of  this  solution, 
and  introduced  a  discontinuity  in  the  saturation-distance  curve  to 
resolve  the  problem.  In  1959,  Fayers  and  Sheldon  (4)  showed  that  the 
introduction  of  the  capillary  term  eliminates  the  multiple-valued 
behaviour  of  the  saturation,  and  in  1976,  Bentsen  (5)  derived  the 
conditions  under  which  the  non-capillary  solution  to  the  immiscible 
displacement  equation  is  an  acceptable  approximation  to  the  transient 
solution. 

Mathematical  models  of  relative  permeabi 1 ity-saturation  and 
capillary  pressure-saturation  curves  are  used  in  numerical  simulators, 
and  must,  of  course,  comply  with  the  experimental  evidence,  as  well 
as,  for  the  capillary  pressure-saturation  relationship  at  least,  with 
some  mathematical  criteria. 

Several  methods  of  capillary  pressure  measurement  have  been 
discussed  in  the  literature.  The  most  commonly  used  are: 

i)  the  diaphragm  method 

ii)  the  mercury  injection  method 

iii)  the  dynamic  method 

iv)  the  centrifugal  technique 
It  is  the  purpose  of  this  study  to  investigate  a  quick 
method,  first  proposed  by  Szabo  (6)  in  1974,  for  determining  the 
capillary  properties  of  consolidated  porous  media.  It  makes  use  of 
a  centrifuge,  and  was  reported  to  be  suitable  for  obtaining  the 
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entire  hysteresis  loop  in  a  matter  of  days.  Drainage  and  imbibition 
cells,  designed  in  principle  after  those  developed  by  Szabo  (6),  were 
constructed  and  used  in  the  experimental  work.  A  new  mathematical 
model  is  proposed  to  fit  the  data  and  is  compared  with  another  model 
previously  studied  by  Bentsen  and  Anli  (7). 


CHAPTER  II 


LITERATURE  REVIEW 


II. 1  CAPILLARY  PRESSURE 

The  capillary  pressure  across  the  interface  between  two 
immiscible  fluids  contained  in  a  capillary  is,  by  definition,  the 
pressure  difference  between  the  non-wetting  and  the  wetting 
phases  (8);  it  can  be  related  to  the  difference  in  level  between  the 
interface  in  the  capillary  and  the  free  fluid-fluid  interface  at 
which  P  =  0  (9),  and  to  the  curvature  of  the  interface  in  the 

V 

capillary  (10,  11).  The  corresponding  expressions  are: 


p 

=  P  -  P 

(2-1) 

c 

nw  w 

p 

c 

=  Ip  ah 

(2-2) 

Pc. 

=  YC 

(2-3) 

where 

P  is  the  pressure  in  the  non-wetting  phase 

P  is  the  pressure  in  the  wetting  phase 

h  is  the  level  of  the  interface  above  or  below  the 

free  fluid-fluid  interface 
Ap"  is  the  average  value  of  the  density  difference 
between  both  fluids  over  the  height  h 
y  is  the  in ter facial  tension 
c  is  the  mean  curvature  of  the  interface 
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These  equations  are  based  on  the  assumptions  that  the  fluids  are 
isothermal  and  in  capillary  equilibrium  (12),  without  any  other 
assumptions  regarding  the  fluids  or  system. 

Because,  in  a  given  pore  and  for  a  given  system,  the  mean 
curvature  is  a  function  of  saturation,  it  may  be  stated  that  a 
functional  relationship  should  exist  between  the  capillary  pressure 
and  the  saturation  (13).  This  will  be  considered  in  a  later 
section. 

It  is  well  known  that  rings  of  isolated  bulk  liquid  are 
formed  during  drainage  and  that  some  of  the  non-wetting  phase  is 
bypassed  during  imbibition.  Hassler  et  al.  (14)  recognized  that 
these  immobile  rings  and  blobs  are  not  in  true  equilibrium  with  the 
main  body  of  fluid  which  surrounds  them,  but  are  in  apparent 
equilibrium;  that  is,  the  changes  within  the  system  are  very  slow. 
Morrow  and  Harris  (15)  have  shown  that  these  isolated  rings  and 
blobs  of  liquid  are  formed  at  the  capillary  pressure  which  existed 
within  the  continuous  bulk  fluid  at  the  time  of  their  formation,  and 
that  they  remain  at  this  pressure  for  a  long  period  of  time,  so  that 
their  pressure  is  considered  to  be  constant.  Second  order  mass 
transfers  such  as  evaporation,  pellicular  fluid  flow,  diffusion,  etc., 
although  they  occur,  are  very  slow  processes  and  are  generally  not 
taken  into  account  in  the  study  of  immiscible  displacement.  Therefore, 
when  bulk  liquid  flow  is  complete,  the  system  is  said  to  be  in 
equilibrium.  It  is  important  to  make  the  distinction  between  true 
and  apparent  equilibrium  and  to  bear  in  mind  the  fact  that  the 
isolated  rings  are  in  apparent  equilibrium;  it  is  neither  possible 


. 
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for  the  desaturation  to  result  from  the  diminution  of  pendular  rings, 
nor  for  these  rings  to  grow  in  the  early  stage  of  imbibition  until 
they  coalesce,  as  it  is  sometimes  suggested  in  the  literature  (16). 

The  fact  that  equilibrium  is  apparent  also  explains  the  existence 
of  an  infinite  number  of  families  of  hysteresis  loops  for  a  given 
system  (17),  and,  if  sufficient  time  were  allowed  for  true 
equilibrium  to  occur,  imbibition  would  cease  only  at  full  saturation 
08). 

White  et  al.(19)  made  an  interesting  study  on  the  physics 
of  desaturation  in  porous  materials,  with  the  following  assumptions: 

i)  gravity  effects  may  be  neglected 

ii)  the  porous  matrix  is  chemically  and  physically  stable 

iii)  the  porous  medium  is  homogeneous 
Figure  II-l  shows  how  they  divided  the  drainage  curve  in  four  sections 
corresponding  to  different  modes  of  desaturation,  for  which  they 
propose  some  mathematical  relationships  between  capillary  pressure 
and  saturation,  based  on  the  assumption  of  relatively  simple  models 
of  porous  media.  These  divisions  are: 

i)  the  boundary  effect  zone,  where  the  capillary  pressure 
is  still  less  than  the  displacement  pressure,  and  for  which  they 
assumed  a  linear  relationship  between  saturation  and  capillary 
pressure. 

ii)  the  primary  transition  zone,  where  desaturation  occurs 
in  a  discrete  manner,  and  for  which  a  logarithmic  relationship  is 
proposed. 


. 


Capillary  pressure  head,  Pc/Y 


Saturation,  S  ,  in  volume  per  volume 


FIGURE  II-l 
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iii)  the  secondary  transition  zone,  where  a  relative 

permeability  to  the  non-wetting  phase  is  established  and  the  effect 

of  pore  isolation  can  be  neglected.  Here  the  authors  assumed, 

without  giving  any  justification,  that  the  frequency  of  pores  that 

can  desaturate  for  an  incremental  change  in  the  length  parameter  R, 

measured  in  terms  of  constriction  widths,  is  a  normal  distribution. 

iv)  the  residual  desaturation  zone,  where  desaturation 

is  due  to  the  retreat  of  the  non-wetting  phase  interface  into  the 

recesses  of  the  pores,  a  mechanism  difficult  to  accept,  because 

pendular  rings  are  not  in  true  equilibrium  with  the  main  body  of 

fluid.  Here  the  relationship  between  S  ,  and  P  is  given  as 

r  w  c  3 


Pd  ? 

‘  Swi^P~) 
c 


(2-4) 


or 


(2-5) 


where 


S  -  S  . 

S  =  w  wi 

1  -  S  . 
wi 

is  the  capillary  displacement  pressure 
5  is  a  constant 

Although  it  is  not  stated  in  their  paper,  £  should  be  greater  than 
one  in  order  for  the  area  under  the  curve  to  be  finite,  a  condition 
which  is  necessary,  because  this  area  represents  the  external  work 
done  on  the  system  during  the  displacement  (20). 
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White  et  al .  (21)  report  some  values  of  £  which  are  less  than  unity, 
and  Eq.(2-4)  is,  therefore,  not  generally  applicable. 

The  discrete  mechanism  of  displacement  taking  place  in  the 
primary  transition  zone  was  first  observed  by  Haines  (22)  and  is 
often  referred  to  as  a  Haines  jump  or  rheon  (23-25).  Such  jumps  are 
observed  when  an  interface  changes  from  an  unstable  to  a  more  stable 
configuration.  The  term  ison  was  adopted (24)  to  define  the  segments 
of  pressure-volume  displacement  curve  between  rheons;  the  isons 
correspond  to  reversible  displacements.  It  is,  then,  evident  that 
such  irreversible  phenomena  as  the  formation  of  pendular  rings  occur 
during  rheons,  as  well  as  the  reunion  of  entrapped  fluid  with  the 
corresponding  bulk  fluid  (26).  It  was  already  mentioned  that  the 
capillary  pressure-saturation  curve  is  time-independent  (apparent 
equilibrium),  and  it  becomes  clear  that  this  is  so  because  rheons 
take  place  on  time  scale  which  is  many  orders  of  magnitude  shorter 
than  that  of  the  macroscopic  displacement  experiment  (27). 

Melrose  (28)  pointed  out  that  the  mechanism  of  non-wetting  phase 
entrapment  during  imbibition  is  different  from  the  trapping  of  the 
wetting  phase  during  drainage;  during  drainage  the  wetting  phase  is 
trapped  in  the  small  pores,  whereas  during  imbibition  the  probability 
for  the  large  pores  to  be  bypassed  is  the  largest.  However,  the 
imbibition  curve  is  also  made  of  a  succession  of  rheons  and  isons, 
and  may  be  considered  as  time-independent  for  practical  purposes. 
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Such  microscopic  studies  of  the  desaturation  and  imbibition 
processes  (29-33),  although  of  interest  for  a  better  understanding  of 
the  physics  involved,  are  primarily  concerned  with  ideal  systems,  and 
have  yet  to  yield  a  sufficent  degree  of  sophistication  so  as  to  provide  a 
capillary  pressure-saturation  relationship  universally  applicable  to  any 
natural  system.  The  literature  is,  therefore,  rich  in  empirical  or 
semi-empirical  relationships,  based  on  experimental  observations, 
which  are  satisfactory  for  particular  systems.  In  1941,  Leverett  (34) 
proposed  his  widely  used  J-function  for  unconsolidated  sands;  it  is 
the  product  of  interfacial  curvature  and  a  property  of  the  rock 
itself,  referred  to  as  "equivalent  circular  diameter  of  the  voids  in 
the  sand"  (35).  It  is  expressed  as 


where 

K  is  the  permeability  of  the  sand 
$  is  the  porosity 


(2-6) 


This  function  was  modified  by  Rose  and  Bruce  (36),  who  introduced 
a  contact  angle  as  follows: 


Brown  (37)  used  this  form  of  the  J-function  to  show  that  the 
correlation  was  improved  by  restricting  its  use  to  data  obtained  from 
a  single  formation,  and  could  be  further  improved  by  dividing 
materials  according  to  their  lithologic  types  in  the  same  formation. 
In  1952,  Burdine  (38)  used  pore-size  distribution  data  to  calculate 


•: 


relative  permeabilities,  and  in  1954,  Corey  (39)  used  this  work  to 

show  that,  for  a  large  number  of  consolidated  porous  rocks,  the 

dimensionless  saturation,  S,  is  a  linear  function  of  1/P  2,  i.e. 

c 


s  =  (/)2  for  Pc  >  Pd 


(2-8) 


where  is  the  displacement  capillary  pressure.  This  relationship 
is,  however,  not  generally  applicable,  because  the  approximations 
on  which  it  is  based  are  valid  only  for  a  particular  pore-size 
distribution  (40).  In  1964,  Brooks  and  Corey  (41),  on  the 
observation  that  log  S  is  a  linear  function  of  log  P  /y  for  numerous 

w 

isotropic  unconsolidated  media,  suggested  the  following  relationship 


S 


(^)x 

c 


for 


(2-9) 


where  A  is  called  the  pore-size  distribution  index  of  the  medium 
and  is  the  negative  of  the  slope  of  the  log-log  plot  mentioned  above 
Very  recently,  this  model  was  used  by  Sigmund  and  McCaffery  (42). 

Any  capillary  pressure-saturation  model  should  explicity 
or  implicitly  allow  for  variations  of  wettability  from  one  system  to 
another.  It  is  therefore  expected  that  A,  in  the  latter  model,  not 
only  reflects  a  characteristic  of  the  porous  medium  itself,  but  also 
varies  according  to  the  pair  of  fluids  used  in  each  experiment 
conducted  on  a  single  porous  medium. 


•  I 


12 


I I. 2  WETTABILITY 

Craig  (43)  defines  wettability  as  "the  tendency  of  one 
fluid  to  spread  on  or  adhere  to  a  solid  in  the  presence  of  other 
immiscible  fluids".  The  behaviours  exhibited  by  various  systems  can 
be  totally  different,  and  much  controversy  still  exists  in  the 
literature  regarding  the  wettability  of  reservoir  rocks.  Some  authors 
(44)  state  that  rocks  are  preferentially  water-wet,  other  (45)  that 
they  are  totally  water-wet,  some  (46)  have  found  certain  cases  of 
preferentially  oil-wet  rocks;  Hassler  et  al .  (47)  showed  the 
wettability  to  be  a  function  of  saturation,  the  rocks  being  water-wet 
at  high  saturation  and  more  oil-wet  at  low  saturation.  They  explain 
this  fact  by  pointing  out  that  the  larger  grains  are  generally  quartz 
and  calcite,  which,  because  of  their  strong  molecular  force  field, 
are  water-wet,  whereas  the  fine  grains  are  generally  clays  and  other 
kinds  of  soft  minerals,  which  do  not  have  a  strong  surface  attraction 
and  exert  significant  capillary  action  upon  the  fluid  boundaries  at 
low  saturation  only.  It  cannot  be  said,  however,  that  soft  minerals 
have  always  a  tendency  to  be  oil-wet;  clays,  in  particular,  are 
water-wet.  The  argument,  therefore,  does  not  seem  to  be  very  strong. 
Kinney  and  Nielsen  (48)  showed  that  the  wettability  of  natural  cores 
depends  on  which  phase  was  first  in  place,  whereas  Kinney  et  al .  (49) 
showed  that  the  wettability  characteristics  were  not  always 
determined  by  the  phase  which  initially  saturated  the  rock.  Finally, 
Alba  (50)  suggested  that  before  oil  migration  occurred  in  the 
reservoirs,  the  rocks  were  probably  water-wet,  but  that  when  oil 
migration  ended,  the  pores  filled  with  oil,  that  is,  the  biggest  pores, 
became  oil-wet.  His  conclusion  is  therefore  contrary  to  that  of 


. 


13 


Hassler  et  al .  (47) . 

Wettability  is  a  very  important  factor  in  oil  production 
and  has  a  major  effect  on  the  characteristics,  such  as  relative 
permeabilities,  of  multiphase  flow  in  porous  media,  due  to  the  very 
large  surface  area  per  unit  volume  of  reservoir  rocks  (51).  The 
relative  mobility  of  oil  and  water  is  partially  a  function  of 
wettability,  and  this  accounts,  among  other  factors,  for  the 
differences  in  relative  permeability-saturation  curves  for  typical 
water-  and  oil-wet  rocks  [Figure  1 1-2] .  For  a  given  oil  viscosity, 
the  mobility  ratio  is  more  favourable  in  the  case  of  a  strongly 
water-wet  rock  than  in  the  case  of  a  strongly  oil-wet  rock,  and  the 
displacement  efficiency  of  the  former  is  higher  than  that  of  the 
latter. 

Wettability  is,  however,  not  easy  to  measure.  It  is 
generally  accepted  that  contact  angles  give  the  true  specification 
of  wettability  (52,  53);  and  some  very  extensive  studies  on  the 
relationship  between  contact  angle  and  capillarity  (54)  and  the  effects 
of  surface  roughness  on  contact  angle  (55)  have  been  reported  in  the 
literature.  These  studies  have  been  made  on  porous  samples  of  PTFE 
(polytetrafluoroethylene)  with  highly  purified  liquids,  and  are 
interesting  as  far  as  the  fundamental  knowledge  and  understanding  of 
phenomena  related  to  contact  angles  are  concerned;  but  it  would  be 
difficult  to  extrapolate  these  results  to  natural  systems,  due  to  the 
heterogeneity  of  mineral  composition  and  the  highly  complicated 
composition  of  crude  oils.  It  has  also  been  recognized  (56-58)  that 
contact  angles  vary  with  time,  due  to  the  adsorption  on  the  solid 
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FIGURE  1 1-2 
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surface  of  polar  components  present  in  minute  quantities  in 
different  crudes.  Contact  angles  also  vary  due  to  deposits  of  asphalt 
or  wax-like  substances  (59)  and  to  other  surface  active  agents  which 
can  also  have  an  effect  on  interfacial  tension.  Hence,  although  still 
widely  used,  a  single  contact  angle  to  characterize  the  wettability 
of  a  particular  system  may  not  always  lead  to  satisfactory  results. 
Thus,  Brown  and  Fatt  (60)  used  the  nuclear  magnetic  relaxation 
method,  instead  of  the  contact  angle,  for  measuring  wettability. 

This  method  is  based  on  the  observation  that  the  adjustment  or 
relaxation  rate  of  fluids  put  in  a  different  magnetic  environment  is 
greatly  affected  by  the  surfaces  of  porous  media.  The  same  year, 

Coley  et  al .  (61)  suggested  that  the  ratio  swj/sor  be  used  as  a 
scale  for  wettability,  because  it  increases  with  wettability.  Other 
tests  have  been  suggested  to  evaluate  wettability.  Amott  (62)  used 
a  method  whereby  he  measured  the  amount  A  of  oil  displaced  by 
imbibition  of  water  in  a  fully  oil-saturated  core,  and  the  amount  B 
of  oil  displaced  by  water  when  centrifuging  the  core  after  free 
imbibition.  Taking  the  ratio  A/ (A  +  B)  gives  an  index  of  wettability: 
if  the  ratio  approaches  one,  the  core  is  strongly  water-wet.  A  similar 
test  is  done  in  which  the  core  is  fully  water-saturated  and  oil  is 
the  displacing  fluid.  If  C  and  D  are  the  amounts  of  water  extracted 
in  the  two  steps  of  the  experiment,  the  ratio  C/(C  +  D)  is  calculated, 
and  a  val  ue  of  one  indicates  that  the  core  is  strongly  oil-wet.  If 
both  ratios  are  equal  to  zero,  there  is  neutral  wettability. 

Another  suggestion  (63)  is  to  evaluate  the  logarithm  of  the 
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ratio  of  the  areas  under  the  drainage  and  imbibition  curves  obtained 
with  the  use  of  a  centrifuge.  The  method  did  not  allow  for  the 
acquisition  of  the  water  imbibition  curve  for  pressures  greater  than 
zero,  nor  the  water  drainage  curve  for  pressures  less  than  zero.  For 
water-wet  systems  the  logarithm  of  the  area  ratio  is  positive,  while 
for  oil -wet  systems  it  is  negative.  This  method  was  reported  by  the 
author  to  be  reliable,  simple  and  fast. 

Wettability  may  have  an  effect  on  the  shape  of  the  capillary 
pressure-saturation  curve.  Fatt  and  Klikoff  (64)  suggested  that 
increasing  the  fractional  oil  wettability  of  a  uniform  grain-size 
sand  pack  lowers  the  capillary  pressure  curve  without  affecting  its 
shape;  but  if  the  grain-size  distribution  is  large,  the  shape  is  also 
affected.  It  may  therefore  be  misleading  to  interpret  capillary 
pressure  curves  in  terms  of  pore  size  distribution,  if  the  wettability 
is  non  uniform.  In  1971,  Morrow  and  Mungan  (65)  reported  some  results, 
obtained  on  a  PTFE  core,  showing  that  the  capillary  pressure  function 
does  not  depend  on  contact  angles  up  to  approximately  49°,  and  that 
it  systematically  decreases  as  contact  angles  increase  from  73°  to 
108°. 

The  actual  knowledge  with  respect  to  wettability  is  still 
sparse  and  does  not  allow  quantitative  analysis  of  the  phenomenon. 

In  particular,  there  is  no  agreement  regarding  the  method  of 
measurement  which  should  be  used.  Contact  angle  as  a  scale  for 
wettability  has  been  widely  challenged  for  the  reasons  referred  to 
above;  another  difficulty  regarding  contact  angle,  as  well  as 
capillary  pressure,  is  the  hysteresis  phenomenon. 
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II. 3  HYSTERESIS 

When  the  state  of  a  physical  system  is  changed  from  A  to  B, 
three  types  of  behaviour  may  be  encountered  (66): 

i)  equilibrium  processes:  if  a  system  is  slowly  changed 
from  state  A  to  state  B,  and  if  the  path  from  B  to  A  is  the  same  at 
any  instant  as  the  path  from  A  to  B,  the  process  is  reversible  or  in 
equilibrium. 

ii)  hysteresis:  if  AB  differs  from  BA,  but  the  paths  are 
stable,  the  system  exhibits  hysteresis. 

iii)  metastability  or  supersaturation:  if  the  paths  are 
not  completely  reproducible,  the  process  is  metastable. 

Hysteresis  is  encountered  in  different  fields  of  science 
and  has  been  extensively  studied  (31,  67-71).  A  very  general  approach 
has  been  taken  by  Everett,  Whitton  and  Smith,  who  postulated  that 
systems  comprising  a  large  number  of  independent  domains,  some  of 
which  might  exhibit  metastability,  exhibit  hysteresis.  These  domains 
are  so  small  that  it  is  not  possible  to  detect  the  behaviour  of  any 
one  of  them,  and  groups  of  similar  domains  are  assumed  to  be  small. 

This  theory  has  been  used  successfully  by  Poulovassil is  (70), 
who  was  able  to  predict  the  path  of  the  primary  drying  scanning  curves 
from  a  knowledge  of  experimental  primary  wetting  scanning  curves; 
these  curves  originate  respectively  from  the  imbibition  and  drainage 
boundary  curves  of  the  hysteresis  loop  and  terminate  at  the  extremes 
of  the  loop.  He  worked  on  constant  structured  discs  and  assumed  that 
advancing  and  receding  contact  angles  were  both  equal  to  zero;  the 
discs  were  made  of  glass  beads  and  were  therefore  strongly  water-wet. 
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Two  basic  assumptions,  however,  limit  the  use  of  this 

theory: 

i)  each  domain  can  be  in  either  of  two  states  only, 
ii)  the  state  of  each  domain  is  independent  of  the  state 
of  its  neighbours;  in  other  words,  no  provision  is  made  for 
cooperative  behaviour. 

Because  the  areas  under  the  capillary  pressure  curves 
correspond  to  the  reversible  work  done  on  or  by  the  system  during 
displacement  (20),  it  is  clear  that  hysteresis  is  a  consequence  of 
the  fact  that  less  energy  is  expended  by  the  system  during  the 
imbibition  process  than  has  been  absorbed  during  drainage.  This 
energy  may  be  stored  by  the  system  in  two  distinct  ways  (72): 

i)  creation  of  an  interface  between  the  continuous 
wetting  phase  and  the  trapped  non-wetting  phase. 

ii)  change  in  wettability  and/or  contact  angle. 

It  should  be  mentioned  here  that  contact  angle  and  capillary  pressure 
hysteresis  are  distinct  phenomena.  The  latter  is  due  to  the 
morphology  of  the  interconnected  pore  space,  within  which  a  definite 
range  of  stable  fluid-fluid  configurations  may  exist.  Both 
interaction  between  pores  (73)  and  a  continuous  range  of  values  for 
the  contact  angles  (74)  have  been  shown  to  exist.  In  fact,  it  is 
the  physics  of  saturation  and  desaturation  itself  which  is 
responsible,  at  least  partially,  for  the  capillary  pressure 
hysteresis:  the  advancing  non-wetting  phase  preferentially 
penetrates  the  large  pores  upon  drainage  and  bypasses  the  small  pores, 
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saturated  with  the  wetting-phase;  whereas  the  advancing  wetting- 
phase  preferentially  penetrates  the  small  pores  upon  imbibition  and 
bypasses  the  larger  ones.  It  is  then  expected  that,  at  any  given 
capillary  pressure,  the  imbibition  curve  will  be  under  the  drainage 
curve,  if  the  wetting-phase  saturation  is  plotted  along  the  abscissa. 

Because  interaction  between  pores  as  well  as  a  finite 
range  of  stable  fluid-fluid  configurations  may  exist,  neither  of  the 
two  basic  assumptions  made  in  the  independent  domain  theory  of 
hysteresis  is  an  acceptable  one,  as  far  as  immiscible  flow  in  porous 
media  is  concerned.  Topp  and  Miller  (75)  showed  that  the  theory  is 
not  generally  applicable. 

II. 4  MEASUREMENT  OF  CAPILLARY  PRESSURE  IN  SMALL  CORE  SAMPLES 

Special  techniques  have  been  devised  in  order  to  measure 
capillary  pressure  in  small  core  samples.  These  are  necessary 
because,  according  to  the  equation 

PC  =  Ap  gh  (2-2) 

high  capillary  pressures,  by  natural  gravitation,  are  obtainable 
only  with  the  use  of  large  cores.  Long  unconsolidated  sand  packs 
are  easy  to  make,  but  long  natural  homogeneous  consolidated  cores 
are  generally  not  available. 

In  1947,  Bruce  and  Welge  (76)  used  the  restored  state 
technique,  which  is  also  known  as  the  displacement  cell  or  diaphragm 
method  (77).  One  of  the  faces  of  the  core  to  be  studied  is 
carefully  machined  flat;  the  core  is  then  cleaned,  saturated  with 
brine  and  placed  in  capillary  contact  with  a  brine-saturated  membrane. 
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whose  average  pore  radius  is  several  orders  of  magnitude  smaller  than 
that  of  the  sample,  so  that  its  displacement  pressure  is  higher  than 
any  pressure  used  in  the  experiment.  The  core  is  placed  on  top  of 
the  disc  and  is  surrounded  by  oil,  whereas  the  other  side  of  the 
diaphragm  is  immersed  in  brine.  The  oil  pressure  is  increased  by 
steps  in  such  a  way  that  time  is  allowed  between  each  pressure 
increment,  so  that  equilibrium  can  be  reached.  When  equilibrium  is 
obtained,  the  saturation  is  measured  and  further  desaturation  is 
carried  on  by  increasing  the  oil  pressure. 

The  same  method  can  be  used,  with  minor  changes,  to  obtain 
the  imbibition  branch  of  the  hysteresis  loop. 

Any  combination  of  fluids  can  be  used  in  this  test,  and 
it  is  generally  considered  to  be  the  best  and  most  reliable  method 
to  obtain  capillary  pressure-saturation  curves.  Its  biggest 
disadvantage  is  the  time  requirement,  which  is  generally  several  weeks 
for  a  single  sample. 

In  1949,  Purcell  (78)  used  the  mercury  injection  method, 
in  which  mercury  is  the  non-wetting  phase  and  mercury  vapour  and  air 
are  the  wetting-phase.  A  clean,  dry  core  is  placed  in  a  chamber  of 
known  volume,  and  a  vacuum  is  drawn  on  the  cell.  The  chamber  is  then 
filled  with  mercury,  and  nitrogen  is  used  to  apply  pressure  at  the 
top  of  the  chamber;  mercury  penetrates  the  core,  and  the  corresponding 
displaced  volume  is  accurately  measured  byrefilling  the  cell.  The 
nitrogen  pressure  is  further  increased,  and  the  incremental  volume 
of  mercury  in  the  core  is  measured.  It  is,  then  possible  to  get  a 
pressure-volume  of  mercury  injected  curve.  Because  of  the  use  of 
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mercury,  the  data  must  be  converted  in  order  to  make  the  comparison 
with  other  systems  possible,  by  taking  into  account  the  differences 
of  interfacial  tensions  and  contact  angles.  The  author  compared  the 
results  obtained  with  mercury  with  those  obtained  on  similar  cores 
with  the  water-air  system.  Assuming  a  capillary  pressure  equation 
of  the  form 

p  =  Zycosjl  (2-10) 

C  K 


a  comparison  of  capillary  pressures  at  corresponding  saturations 
leads  to  the  ratio: 


P 


c  Hg 


(ycose) 


c  air/water 


(ycose ) 


H£ 


air/water 


(2-11) 


The  following  values  were  assumed  to  be  acceptable: 

yu  =  480  mN/m 

Hg 

y  ,  =70  mN/m 

ywater 

e u  /  i  •  i  =  140° 

Hg/sol id 

0 water/sol  id  ^ 

and  consequently, 


c  Hg 


c  water/air 


=  5.25  -  5. 


Fairly  good  agreement  was  obtained  for  both  systems.  There  are  some 
doubts,  however,  concerning  the  validity  of  the  assumptions  made, 
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and  in  1951,  Brown  (37)  found  that  by  using  a  conversion  factor  in  the 
range  5.4  to  8.3  for  different  cores,  he  could  get  almost  perfect 
agreement  with  other  data. 

The  main  advantages  of  the  mercury  injection  method  are 
that  it  is  very  fast,  it  allows  very  high  pressures  to  be  obtained, 
and  any  kind  of  irregularly  shaped  sample  may  be  used.  However,  the 
interpretation  of  data  is  uncertain,  and  the  test  results  in  the 
permanent  loss  of  the  samples. 

In  1951,  Brown  (79)  pointed  out  that,  in  spite  of  the  fact 
that  both  the  diaphragm  and  the  mercury  injection  methods  give  the 
capillary  pressure  curves  measured  under  static  equilibrium,  the 
results  thus  obtained  were  used  to  solve  dynamic  flow  problems.  He 
modified,  therefore,  a  test  previously  devised  by  Hassler  (80)  for 
permeability  measurements,  and  rendered  it  suitable  for  capillary 
pressure  measurements  under  dynamic  equilibrium.  In  this  test,  a 
simultaneous  steady  flow  of  oil  and  gas  is  allowed  through  one  end 
of  the  core.  The  pressure  drop  in  the  oil  phase  is  measured  through 
specially  designed  semi -permeable  membranes  placed  at  the  inlet  and 
outlet  ends  of  the  core,  and  the  flow  of  gas  is  monitored  so  that  the 
pressure  difference  in  both  phases  -  that  is,  the  capillary  pressure  - 
is  the  same  on  both  sides.  As  a  constant  capillary  pressure  is 
established  along  the  core,  a  constant  saturation  also  exists.  By 
varying  the  flow  rates,  a  complete  curve  may  be  obtained.  The 
agreement  between  dynamic  and  static  methods  was  found  to  be  very  good; 

Brown  (37)  concluded,  therefore,  that  it  is  justified  to  use  data 
obtained  under  static  equilibrium  for  dynamic  flow  problems.  Such 


. 
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a  conclusion,  however,  is  not  totally  justified;  it  is  true  that  a 
dynamic  method  is  closer  to  actual  flow  and  displacement  systems 
than  a  static  method,  and  automatically  corrects,  at  least  partially, 
for  such  related  factors  as  the  difference  between  static  and 
dynamic  interfacial  tensions  or  receding  and  advancing  contact  angles. 
However,  it  should  be  borne  in  mind  that  this  method  is  also  a  method 
in  which  equilibrium  is  established  before  measurements  can  be  made. 

An  interesting  feature  of  this  test  is  that  the  measured  saturations 
are  true  and  not  average  saturations,  provided  the  core  is  homogeneous. 
It  is,  however,  slow  and  requires  much  care. 

It  was  deemed  necessary  to  devise  a  method  whereby 
equilibrium  would  obtain  quickly.  The  centrifuge  method  was  probably 
used  for  the  first  time  by  McCullough  et  al . (81 ) ,  but  the  first 
extensive  consideration  of  its  mathematical  aspect  is  due  to  Hassler 
and  Brunner  (82),  who  also  used  an  entirely  different  core-holding 
cell.  Referring  to  Eq . (2-2 ) ,  it  is  obvious  that  high  capillary 
pressures  can  be  obtained  either  with  very  long  core  samples  in  the 
earth  field  of  gravity,  or  with  an  increased  field  of  gravity  on 
small  core  samples.  The  centrifuge  method  makes  use  of  the  latter 
option. 

The  capillary  pressure  may  be  determined  provided  the  speed 
of  rotation,  the  density  difference  between  the  fluids,  and  a  few 
geometrical  parameters,  resulting  from  the  design  of  the  equipment 
used,  are  known;  the  average  saturation  within  the  core  is  derived 
from  a  material  balance.  It  is  thus  possible  to  get  a  curve  relating 
the  capillary  pressure  at  any  point  in  the  core  to  the  average 
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saturation  within  the  core  at  this  capillary  pressure.  In  the 
following  derivation,  two  assumptions  are  made: 


i)  the  capillary  pressure  is  zero  at  the  outer  end  of 
the  core. 

ii)  the  variation  of  acceleration  along  the  core  is 
negl igible. 


Let 


L  be  the  length  of  the  core 


Ap  the  density  difference  between  the  fluids  used 
a  the  acceleration 

h  the  distance  above  the  outer  end  of  the  core  at  which 
the  capillary  pressure  is  calculated. 

The  saturation  Swalong  the  core  is  a  function  of  the 


capillary  pressure,  and  the  average  saturation  is  given  by 


(2-12) 


0 

By  choosing  Pc1-=  ApgL  as  the  variable,  Eq .  (2-12)  may  be  written  as 


P  . 


(2-13) 


0 


from  which  it  follows  that 


S  (P  •)  =  -tJ-  (S  P  •) 

wv  ci '  dP  .  K  vj  ci ' 
ci 


(2-14) 


Hassler  and  Brunner  (83)  used  graphical  differentiation  to 


obtain  the  capillary  pressure-saturation  curve. 
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Let  r.  and  rg  be,  respectively,  the  distances  of  the  inner 
and  outer  ends  of  the  core  to  the  axis  of  rotation,  and  oo  the  angular 
velocity  of  the  centrifuge.  If  the  variation  of  the  centrifugal  field 
is  taken  into  account,  the  following  equations  apply: 


Pc  =  Apw2  (rg  -  2")  h 


Pci  =  1  Apco2  (r2  -  r?)  =  Apw2  (r0  -  j)  L 


dP. 


dh  = 


wVi  -  P77  (re  - 

Therefore,  Eq . (2-12)  may  be  written  as  follows: 

P 


or 


°w 


ci  Sw(Pc)  dP( 


ApLuT 


Vre  -  F-  (i  -  '?> 


Cl 


Pci 


c--  'e  2 

w  ~77’^7 


ypc)  dp< 


0 


1  Cl 


r? 

(1  -  -) 


Therefore, 


P  -S,, 
ci  w 


-  \  o  +  ?■) 


Cl 


ypc) dp, 


0  Vi  -  pS1  -  ^ 

ci  r0 


and 


P  . 
ci 


P  .S 
ci  w 


=  cos 


.2  L 


Sw(Pc>  dPc 


V1  - 


sin2f 


ci 


ri/re' 


(2-15) 

(2-16) 

(2-17) 


(2-18) 


(2-19) 


(2-20) 


(2-21) 


where  cos  f 
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Equation  (2-21)  cannot  be  solved  simply  for  the  function 
S  ,  and  a  method  of  successive  approximations  was  proposed  by  the 
authors.  This  method,  however,  is  very  tedious,  and  they  reported 
that  for  values  of  the  ratio  r./r  as  small  as  0.7,  the  approximation 
given  by  Eq . (2-14)  is  usually  sufficient. 

Few  data  were  given  by  Hassler  and  Brunner  (82).  In  1951, 
Slobod  et  al .  (84)  reported  results  which  showed  the  method  to  be 
rapid  and  reliable;  high  capillary  pressures  were  obtained  as  well  as 
values  of  connate  water  and  residual  oil  saturations. 

At  this  point  in  time,  however,  only  one  branch  of  the 
hysteresis  loop,  namely  drainage,  could  be  obtained  by  this  method, 
and  according  to  Donaldson  et  al.  (63): 

"the  water  imbibition  curve  for  pressures  greater 
than  zero  and  the  water  drainage  curve  for  pressures 
less  than  zero  cannot  be  obtained  by  the  centrifugal 
technique" . 

In  1974,  Szabo  (6)  devised  three  different  core-holding 
cells,  whereby  it  is  possible  to  measure  imbibition  capillary 
pressure  using  the  centrifugal  technique.  All  of  them  meet  the 
following  necessary  requirements: 

i)  the  direction  of  phase  exchange  occurs  as  a  result  of 
the  natural  fluid  differences;  that  is,  the  water  enters  the  sample 
moving  off  the  rotary  axis. 

ii)  the  quantity  of  imbibed  water  is  controllable. 

To  achieve  this,  a  semi-permeabl e  filter  disc,  which  separates  the 
water  reservoir  from  the  sample,  is  placed  in  capillary  contact  with 
the  inner  end  of  the  sample.  The  water  reservoir  is  above  the  filter 
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disc  and  is  opened  to  the  atmosphere.  The  disc  is  water  saturated, 
and  the  sample  is  surrounded  by  the  oil  phase.  When  starting  the 
imbibition,  a  speed  is  chosen  that  will  allow  the  capillary  pressure 
inside  the  core  to  be  somewhat  smaller  than  the  capillary  pressure 
created  at  the  end  of  the  drainage  experiment.  When  equilibrium  is 
obtained,  the  average  saturation  is  determined  by  weighing  or  by 
material  balance,  and  a  lower  speed  is  selected  to  allow  further 
imbibition.  A  decrease  in  speed  results  in  a  decrease  in  capillary 
pressure,  and  therefore  imbibition  occurs.  The  time  necessary  to 
stop  the  centrifuge  is  so  small  as  compared  to  the  time  necessary  for 
the  fluids  inside  the  core  to  reach  equilibrium,  that  Szabo  (85)  assumed 
it  does  not  allow  further  imbibition  within  the  core  to  occur. 

The  first  technique  consists  of  measuring  the  capillary 
pressure  within  a  small  core  sample.  Imbibition  as  well  as  drainage 
could  be  obtained,  but  the  design  of  the  cell  is  such  that  both  oil 
and  water  columns  are  large  and  limit  the  maximum  speed  of  rotation 
required.  Moreover,  the  cell  must  be  disassembled  after  each  run, 
and  the  method  is  therefore  slow  and  tedious. 

A  second  method  was  designed  to  measure  the  imbibition 
capillary  pressure  in  cores  of  intermediate  length.  Because  of  the 
length  of  the  samples  used  in  this  method,  the  difference  in  the 
capillary  pressures  at  both  ends  of  the  samples  is  significant,  and 
the  saturation  varies  along  their  length.  An  average  saturation  is 
thus  obtained,  and  the  author  proposed  the  use  of  graphical 
differentiation  in  order  to  calculate  the  saturation  at  the  inner 
end  of  the  core.  Provided  an  acceptable  capillary  pressure  model  is 
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found,  a  numerical  solution  to  this  problem  would  be  preferable  (86). 

The  third  method  suitable  for  drainage  as  well  as 
imbibition  capillary  pressure  measurements  was  referred  to  as  the 
"quick  method"  by  the  author.  It  consists  of  using  a  long  core 
sample  cut  into  slices,  which  are  held  in  capillary  contact  by  discs 
of  filter  paper  placed  in  between  each  slice.  Upon  reaching 
equilibrium  the  average  saturation  within  each  slice  is  measured  by 
weighing,  and  Szabo  (87)  plotted  the  capillary  pressure  calculated 
at  the  center  of  each  slice  against  the  average  saturation  of  the 
corresponding  slice. 


CHAPTER  III 


STATEMENT  OF  PROBLEM 


Capillary  pressure  is  a  very  important  phenomenon,  and  it 
is  deemed  necessary  to  be  able  to  obtain  reliable  capillary  pressure- 
saturation  data  in  short  periods  of  time.  The  centrifuge  method 
seems  to  be  the  only  one  which  satisfies  this  time  requirement.  This 
method,  however,  is  not  direct,  and  it  is  necessary  to  convert  the 
experimental  data  into  a  true  capillary  pressure  curve.  It  is, 
therefore,  important  to  get  a  good  understanding  of  what  the 
centrifuge  data  really  represent. 

The  quick  method  proposed  by  Szabo  (88)  was  selected  to 
obtain  capillary  pressure-saturation  curves.  Reasons  for  the  choice 
of  this  technique  were: 

i)  it  seems  suitable  to  obtain  drainage  and  imbibition 
curves. 

ii)  several  points  are  obtained  at  once,  and  therefore 
the  method  is  fast. 

iii)  a  saturation  close  to  the  irreducible  saturation  may 
be  obtained  more  easily  than  if  a  core  of 
intermediate  length  were  used, 

iv)  if  the  slices  are  small,  there  is  a  lesser  effect  of 
gravity  as  compared  to  conventional  methods; 
therefore,  the  experimental  curve  should  be  closer 
to  the  true  curve. 
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v)  the  end  effect  is  expected  to  be  of  lesser  importance, 
as  compared  to  that  obtained  in  the  conventional 
centrifugal  technique,  because  the  range  of  rotational 
speeds  for  the  quick  method  is  smaller  than  that 
necessary  to  obtain  the  capillary  pressure  curve  with 
one  core  of  small  or  intermediate  length. 

Most  capillary  pressure  studies  are  restrained  to  the 
drainage  process,  whereas  waterflooding  is  mainly  concerned  with  the 
imbibition  mechanism.  Attempts  were  therefore  made  to  obtain 
imbibition  data  using  the  centrifugal  method. 

Finally,  a  new  mathematical  model  relating  capillary  pressure 
and  saturation  was  studied  to  convert  the  experimental  data  into  a  true 
capillary  pressure  curve;  this  was  rendered  necessary  due  to  the  fact 
that  the  average  saturation  in  the  core,  and  not  the  saturation 
corresponding  to  the  calculated  capillary  pressure,  was  measured. 


CHAPTER  IV 


CAPILLARY  PRESSURE  MODEL 


IV. 1  DRAINAGE 

IV. 1.1  Drainage  Model 

Because  the  capillary  pressure  at  any  point  is  directly 
related  to  the  mean  curvature  of  the  interface  at  that  point,  and 
because,  for  a  given  system,  the  mean  curvature  is  a  function  of  the 
saturation  (9),  it  should  be  possible  to  express  the  capillary 
pressure  as  a  function  of  saturation. 

Two  kinds  of  drainage  capillary  pressure  curves  are  found 
in  the  literature  [Figure  I V- 1 ] .  In  his  study  of  the  role  of 

gravity  in  capillary  pressure  measurements,  Szabo  (89)  showed  that 
the  concavity  of  curve  2,  near  the  maximum  saturation,  is  due  to  the 
effect  of  gravity  on  cores  of  finite  dimensions  and  is  not  a 
characteristic  of  the  drainage  capillary  pressure-saturation 
relationship  itself.  Curve  1,  on  the  contrary,  represents  the 
capillary  pressure-saturation  curve  for  a  very  thin  porous  disc;  it 
approaches,  therefore,  the  true  capillary  pressure  curve,  because  it 
relates,  in  the  limit,  the  capillary  pressure  and  the  saturation  in 
the  corresponding  plane,  excluding  any  gravity  effect.  Therefore, 
a  mathematical  model  plotting  as  curve  1  is  thought  to  be  suitable. 

Several  empirical  or  semi-empirical  mathematical  models 
have  been  used  in  the  literature  (7,  19,  41,  90,  91);  some  of 
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Capillary  pressure 
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Water  saturation 


FIGURE  IV-1 


CAPILLARY  PRESSURE  CURVES  ENCOUNTERED 
IN  THE  LITERATURE  (89) 
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them  (90)  are,  a  priori,  not  acceptable,  because  they  fail  to  meet 

several  basic  requirements.  These  requirements  are  the  consequence 

of  the  physics  involved  as  well  as  of  the  necessary  compatibility  of 

the  models  with  the  displacement  equation  of  one  fluid  by  another, 

and  may  be  summarized  as  follows  (7): 

i)  the  slope  of  the  drainage  capillary  pressure  curve 

must  tend  to  minus  infinity  when  the  water  saturation  approaches  the 

irreducible  water  saturation. 

ii)  the  drainage  capillary  pressure  curve  must  tend  to  a 

fixed,  finite  value,  as  the  water  saturation  approaches  one. 

iii)  the  area  under  the  drainage  capillary  pressure- 

saturation  curve  must  finite,  as  this  area  corresponds  to  the 

reversible  external  work.  Moreover,  any  functional  relationship 

postulated  must  of  course,  agree  with  experimental  evidence. 

The  first  of  these  modelling  conditions  was  proposed  by 

Fayers  and  Sheldon  (92),  who  assumed  that  dP  /dS  behaves  like  -1/S 

>■> 

near  the  irreducible  water  saturation,  where  S  is  the  saturation 

S  -  S  . 

in  dimensionless  form,  S  =  ^  _  SW1-  .  On  this  basis,  Anli  (93) 

\i 

derived  a  model  defined  in  the  domain 

0  <  S  <  1 
Pd  <  P  <  co 

where  Pd  is  the  displacement  pressure  and  has  a  finite  value 
according  to  the  second  requirement.  Assuming 

dP 

c  a 

dS  "  S 


(4-1) 
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and  integrating  with  respect  to  S  in  the  domain  previously  defined, 
yields 


P_  -  -a  in  S  +  P  , 
c  d 


where 


a  is  a  parameter  which  includes  the  effects  of 


wettability,  interfacial  tension  and  pore-size 
distribution,  and  has  the  dimensions  of  pressure. 
This  model  can  be  written  in  dimensionless  form  as  follows: 


(4-2) 


it  =  -J,n  S 


(4-3) 


where 


TT 


P  - 
c  d 


is  a  dimensionless  capillary  pressure. 


Although  it  meets  all  of  the  criteria  mentioned  above, 
this  model  suffers  from  the  fact  that,  as  the  water  saturation 
approaches  the  irreducible  saturation,  the  capillary  pressure  goes 
to  plus  infinity.  To  remedy  this  situation,  Eg. (4-1)  was  modified 


to 


dP( 

dS 


O 

~n 


n  i  1 


(4-4) 


which,  upon  integration,  yields 


a  (1  -  S1_n)  +  P 


1-n 


(4-5) 


Equation  (4-5)  may  be  written  in  dimensionless  form  as 

.1-n 


7T 


1  -  s 


(4-6) 
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where 


TT 


P  -  P  . 
c  d 


is  a  dimensionless  capillary  pressure. 


This  modification  does  not,  as  it  seems  to,  violate  the  condition 
imposed  by  Fayers  and  Sheldon  (92),  and  is  therefore  valid;  this  is 
shown  in  Appendix  A. 

It  is  immediately  evident  that  this  model  meets  the  first 
and  second  requirements.  The  area  under  the  curve  is  given  by 


area  = 


1  -n 


S  - 


$2-n 

2-n 


l 

+  PH 
0  d 


(4-7) 


and,  therefore,  n  must  be  strictly  less  than  2  in  order  for  the 
third  requirement  to  be  met.  However,  for  1  <  n  <  2,  the  capillary 
pressure  goes  to  plus  infinity  as  the  water  saturation  approaches  the 
irreducible  saturation;  therefore,  the  condition  that  0  <  n  <  1,  was 
imposed  on  the  model.  The  agreement  with  the  experimental  evidence 
is  one  of  the  objects  of  the  experimental  study  and  will  be 
considered  later  on. 

Equation  (4-5)  may  also  be  written  as 


1-n 


wi 


+  d  -  swi>  (p- 


a 


1-n 


P  +  v 
c  1  -n 


(4-8) 


The  experimentally  measured  water  saturation  is,  in  fact, 
an  average  saturation,  because  a  core  has  a  definite  thickness, 
whereas  a  relationship  such  as  expressed  by  Eq . (4-8)  relates  the 
capillary  pressure  in  a  certain  plane  in  the  core  to  the  saturation 
in  the  same  plane;  this  equation  must,  therefore,  be  adapted  to  the 
kind  of  experimental  data  obtained. 


. 
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IV. 1.2  Application  of  the  Model  to  the  Drainage  Data  Obtained  with 
the  Centrifuge  and  the  Multi-Slice  Core  Technique 

The  two  directly  available  data,  with  this  technique,  are 
the  average  saturation  in  each  slice  and  the  capillary  pressure  at 
any  point  along  the  core.  Assuming  that  the  capillary  pressure  is 
zero  at  the  outer  end  of  the  sample  (83),  the  capillary  pressure  is 
given  by 

Pc  =  Aapu2  (re  -  |)  h  (4-9) 


Let  us  consider  the  core  to  be  made  of  a  superposition 
of  q  slices  in  capillary  contact,  and  let  us  number  the  contact 

t  h 

planes  as  shown  in  Figure  IV-2.  The  average  saturation  in  the  j — 
slice  is 


(4-10) 


where  H.  and  h.  are  defined  in  Figure  IV-2. 

3  J 

Combining  Eq.(4-8),  Eq.(4-9)  and  Eq . ( 4-1 0 )  yields 


(4-11) 


where  A  is  the  unit  conversion  factor. 

Equation  (4-11)  is  valid  only  for  the  slices  for  which  >  h*, 
where  h*  is  the  location  of  the  plane  at  which  Pc  =  Pd-  If  h*  is 
located  in  the  g—  slice,  the  average  saturation,  S^,  of  this 


slice  is  given  by 


I ' 
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Multi -si ice  core 


FIGURE  IV-2 
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The  model  proposed  by  Bentsen  and  Anli  (7)  was  also 
used  for  comparison  purposes.  In  this  case,  Eq.(4-2)  must  be 
rewritten  as 


+  (1  -  S 


(4-12) 


(4-13) 


(4-14) 


(4-15) 


Equations  (4-11)  and  (4-12)  are  derived  in  Appendix  B.  These 

equations  cannot  be  solved  analytically  and  were  solved  numerically. 

It  is  evident  that  two  kinds  of  experimental  curves  may  be 

obtained.  On  the  one  hand,  one  may  plot  P.  against  S  . ;  on  the  other 

J  wj 

hand,  one  may  plot  P.  against  S  • ,  where 

J  wj 


(4-16) 


If  the  logarithmic  model  is  used,  then 


S  .  can  be  obtained  from  the  experimental  results  by  Eq . ( 4-1 9 ) : 
wj 
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where 

V- j  is  the  volume  of  water  in  the  i—  slice  at  the  end  of 
the  experiment. 

Vpt  is  the  pore  volume  of  the  i—  slice. 


IV. 2  IMBIBITION 

IV. 2.1  Imbibition  Model 

Assuming  the  imbibition  curve  to  have  the  shape  represented 
in  Figure  IV-3,  the  two  previously  defined  models  were  supposed  to 
be  suitable  to  represent  the  imbibition  capillary  pressure-saturation 
relationship  with,  however,  the  following  modifications: 


S  = 


S  -  S  . 
w  wi 


1  -  S  .  -  S 
wi  or 


(4-20) 


and 


Pd  -  0 


Hence,  Eq . ( 4-2 )  and  Eq . (4-5)  become,  respectively, 


P  =  -  a  £n  S 

c 


(4-21) 


and 


a  /-,  rl-n 


c  1-n 


(1  -  S'"'1) 


(4-22) 


IV. 2. 2  Application  of  the  Model  to  the  Imbibition  Data  Obtained  with 
the  Centrifuge  and  the  Multi-Slice  Core  Technique 

The  experimental  capillary  pressure  is  calculated  by  the 
following  expression  (87),  derived  in  Appendix  D: 


% 


Capillary  pressure 


FIGURE  IV-3 


IMBIBITION  CAPILLARY  PRESSURE  CURVE 
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Aw2 


-[(re-h)2(Po-Pw)  +  r2Pw-r2po]  (4-23) 


Combining  Eg.  (4-10),  Eq.( 4-20) ,  Eq.  (4-21 )  and  Eq .  (4-23)  yields 


s  .  -  S  .  + 

wj  wi 


1  -  s  .  -  s 

_ wi  or 

H  •Ato2/Ap' 


exp  (-  -±) 


W 


dPc  (4-24) 


Am2 


—  +  r2p  -  r2p 
w  w  oMo 


or. 


S  •  -  S  .  + 

wj  wi 


dPc  (4-25) 


where  PQ  is  the  capillary  pressure  at  the  outer  end  of  the  core. 
If  Eq . ( 4-22 )  is  used  instead  of  Eq . (4-21 ) ,  then 


S  • 
wj 


S  .  + 
wi 


H  .Aoo2/Ap" 


dP 


+  r2p  -  r2p 
Aw2  w  w  00 


(4-26) 


or, 


1  -  S  .  -  S  i 

r=  =  S  .  +  - 21 - 2H  (I=H) 

WJ  W1  h.AcoW  0 


(_£_  .  p  ) 

M-n  V 


1 

T=n 


~2P 

— —  +  r2p  -  r2p 
Aw2  w  w  0  0 


dP 


(4-27) 


Equation  (4-24),  Eq . (4-25) ,  Eq . (4-26)  and  Eq . (4-27)  are  derived 
in  Appendix  C. 
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CHAPTER  V 


EXPERIMENTAL  EQUIPMENT  AND  PROCEDURE 

The  experimental  procedure  used  in  this  study  is  similar 
to  that  used  by  Szabo  (88);  he  referred  to  this  method  as  the  "quick 
method  for  capillary  pressure  measurement".  A  long  core  sample,  cut 
into  slices,  is  rotated  in  a  centrifuge.  The  slices,  which  are 
supposed  to  be  homogeneous  and  to  have  the  same  permeabilities,  are 
maintained  in  capillary  contact  with  the  aid  of  Kleenex  paper  placed 
in  between  each  slice;  moreover,  a  coil  spring  acts  on  top  of  the 
core  holder  which  forces  the  slices  against  each  other.  Several 
authors  (88,  94,  95)  have  used  this  procedure  to  maintain  samples  in 
capillary  contact;  it  does  not  appear  to  affect  the  capillary 
pressure  distribution,  and  consequently,  does  not  appear  to  affect 
the  saturation  distribution  within  the  core. 

The  advantage  of  the  quick  method  for  capillary  pressure 
measurement,  as  compared  to  the  conventional  use  of  the  centrifuge 
in  the  study  of  capillary  pressure,  lies  in  the  reduction  of  time 
which  is  required  to  obtain  the  whole  drainage  capillary  pressure 
curve.  Moreover,  Szabo  (85)  presents  some  imbibition  curves  obtained 
with  the  same  method,  and  according  to  his  work,  it  seems  possible 
to  obtain  the  whole  hysteresis  loop  in  a  matter  of  days. 
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V.l  THE  CENTRIFUGE 

An  IEC  PR-6000  centrifuge  was  used  in  this  study,  together 
with  a  model  276  four-place  rotor  for  which  the  maximum  authorized 
rotational  speed  is  3000  rpm.  This  centrifuge  is  provided  with  a 
feedback  speed  control,  whereby  the  selected  speed  is  kept  constant 
within  ±  1.5%,  and  a  temperature  control  device,  for  temperatures 
ranging  from  -10°C  to  +30°C.  It  also  features  a  flexible  drive 
shaft  which  eliminates  the  need  for  critical  balancing  of  the  cells. 

A  window  was  cut  in  the  cover  of  the  centrifuge  to  allow  accurate 
measurement  of  the  speed  of  rotation  with  a  stroboscope,  because  the 
scale  of  the  electric  tachometer  provided  on  the  centrifuge  was 
graduated  in  100  rpm  increments,  and  was  too  coarse  for  the  purpose 
of  this  study. 

V. 2  DRAINAGE 

V.2.1  Drainage  Capillary  Pressure  Cell 

The  drainage  capillary  pressure  cell  was  designed  to 
receive  large  cylindrical  cores,  about  7.2  cm  long  and  2.5  cm  in 
diameter.  A  cross-section  of  the  cell  in  actual  dimensions  is  given 
in  Figure  V-l.  The  cell  body  and  the  core  holder  were  made  of  lucite 
to  allow  for  visual  examination;  this,  however,  has  the  disadvantage 
of  limiting  seriously  the  choice  of  fluids  which  can  be  used  in  the 
experiments.  The  core  holder  was  grooved  vertically  and  bored  with 
a  series  of  holes  at  its  top  to  prevent  the  presence  of  air  bubbles 
around  the  core.  The  top  and  bottom  caps,  as  well  as  the  core 
support  disc,  were  made  of  aluminum;  the  core  support  disc  was 


. 


FIGURE  V-l 


DRAINAGE  CAPILLARY  PRESSURE  CELL 
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radially  and  circularly  grooved  to  allow  free  circulation  of  fluids 
from  the  core  into  the  pipette,  and  from  the  pipette  into  the  cell 
body.  The  pipette,  provided  to  receive  the  water  expelled  from  the 
core,  was  attached  to  the  core  support  disc  by  means  of  a  seal  o-ring 
and  was  made  of  glass  able  to  withstand  an  inside  excess  pressure  of 
16  MPa.  The  bottom  plug  of  the  pipette  was  made  of  teflon  and  sealed 
by  means  of  a  seal  o-ring.  The  arm  and  tube  bodies  were  made  of  mild 
steel  for  a  proper  location  of  the  center  of  gravity.  Two  such 
equally  weighed  cells  were  constructed. 

V. 2.2  Core  and  Cell  Preparation 

The  slices  of  the  core  under  study  were  first  cleaned  with 
toluene  in  a  soxhlet  extractor,  dried  at  150°C,  and  heated  at  500°C 
for  about  10  to  12  hours  to  remove  any  organic  residuals  from  the  pore 
space.  The  dry  slices  were  then  weighed,  placed  under  vacuum,  and 
saturated  with  de-aerated  distilled  water.  The  saturated  slices  were 
weighed  in  water  to  prevent  evaporation,  and  the  mass  of  water  present 
in  each  slice  was  calculated.  It  should  be  pointed  out  that  this 
method  might  lead  to  erroneous  values  of  water  saturation  in  the  case 
where  cores  containing  large  quantities  of  water  of  crystallization 
are  studied. 

Because  each  slice  was  slightly  different  from  the  others, 
it  was  very  important  to  characterize  them  and  to  distinguish 
between  the  inlet  and  outlet  faces.  They  were  therefore  numbered, 
so  that  no  confusion  could  be  possible. 

After  assembling  the  drainage  capillary  pressure  cells, 
the  pipettes  were  filled  with  oil.  Excess  water  was  blotted  from  the 


. 
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surface  of  each  slice,  and  the  cores  were  made  by  alternating  the 
slices,  according  to  their  numbers,  and  small  discs  of  Kleenex  paper 
having  the  diameter  of  the  samples.  The  cores  were  then  inserted  in 
the  core  holders  and  placed  in  the  cells;  these  were  then  filled  with 
oil,  so  that  the  cores  were  totally  immersed. 

A  speed  of  rotation  was  selected,  and  the  drainage 
experiment  was  started.  When  apparent  equilibrium  was  reached,  that 
is,  when  no  more  water  was  expelled  from  the  core  over  a  sufficient 
period  of  time  (10  to  12  hours),  the  centrifuge  was  stopped  and  the 
slices  were  weighed  in  oil  to  prevent  evaporation. 

V.2.3  Drainage  Capillary  Pressure  and  Saturation  Measurement 

The  speed  of  rotation  was  accurately  measured  with  a 
stroboscope,  and  the  densities  of  the  fluids  were  determined  with  an 
AP  Paar  DMA60-Digital  Density  Meter.  The  capillary  pressure  is  given 
by 

Pc  =  1.097  x  10~6  Aprn2  (re  -  h  (5-1) 

where 

pc>  kPa 
Ap,  g/cc 
a  ,  rpm 

r  is  the  distance  from  the  axis  of  rotation  to  the  bottom 
e 

of  the  core,  cm 

h  is  the  distance  from  the  bottom  of  the  core  to  the 
point  where  P  is  calculated,  cm. 

This  equation  is  based  on  the  assumption  that  the  capillary  pressure 
vanishes  at  the  outlet  end  of  the  core. 
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The  average  saturation  in  each  slice  is 


where 


(5-2) 


is  the  mass  of  water  in  the  100%-water-saturated 
slice,  g 

is  the  mass  of  water  and  oil  contained  in  the  slice 
at  the  end  of  the  run,  g 
Vp  is  the  pore  volume  of  the  slice,  cm3 


Equations 


(5-1)  and  (5-2)  are  derived  in  Appendix  D. 


V. 3  IMBIBITION 

V.3.1  Imbibition  Capillary  Pressure  Cell 

The  drainage  capillary  pressure  cell  was  designed  so  that 
it  could  be  used  to  study  imbibition,  provided  the  core  holder  were 
changed.  As  already  mentioned,  it  is  necessary  for  the  water  to  move 
off  the  rotary  axis  in  order  to  take  advantage  of  the  increased 
gravitational  field;  thus,  the  outlet  end  of  the  core  should  not  be 
in  contact  with  the  water  reservoir.  Nor  can  a  situation  in  which  the 
core  is  immersed  in  water  be  considered,  because,  in  this  case,  the 
amount  of  imbibition  could  not  be  regulated.  Therefore,  the  water 
reservoir  must  be  placed  at  the  top  of  the  core.  Because  it  is 
desired  to  study  the  capillary  pressure  curve  in  the  positive 
capillary  pressure  domain,  the  system  must  be  such  that,  at  all  points, 
the  pressure  in  the  oil  phase  is  higher  than  or  equal  to  the  pressure 


49 


in  the  water  phase.  This  is  true,  in  particular,  at  the  top  of  the 
core,  and  therefore  a  semi -permeable  membrane  must  be  used  to  isolate 
the  water  reservoir  from  the  oil  reservoir;  this  membrane  must  be  water- 
wet  and  such  that  its  displacement  pressure  is  higher  than  any 
capillary  pressure  reached  in  the  course  of  the  experiment.  The 
imbibition  capillary  pressure  cell  is  shown  in  actual  dimensions  in 
Figure  V- 2.  The  water  reservoir  was  made  of  lucite  to  allow  for 
visual  examination  of  the  water  level,  and  the  core  holder  was  made 
of  aluminum. 

Many  problems  were  encountered  in  preventing  leakage  of  oil 
into  the  water  reservoir.  Because  it  was  necessary  for  the  membrane 
to  be  permeable  to  water  only,  it  was,  indeed,  deemed  important  to  be 
able  to  clean  it  and  to  heat  it  at  high  temperature,  in  order  to 
eliminate  any  organic  material  from  the  pore  space  and  render  it 
strongly  water-wet.  Different  solutions  have  been  envisaged,  among 
which  a  tight  fit  seemed  to  be  the  best.  Lucite  and  titanium  rings 
were  precision-machined  to  shrink-fit  the  membrane,  but,  although 
much  care  was  taken  to  machine  the  rings  as  well  as  the  membranes,  no 
fit  was  able  to  withstand  an  excess  pressure  of  more  than  28  kPa.  A 
second  possible  solution  was  to  cast  a  low  melting-point  metal  around 
the  membrane;  many  attempts  were  made  with  aluminum,  but  the  presence 
of  air  in  the  pore  space  of  the  disc  always  resulted  in  the  formation 
of  a  small  air  bubble  at  some  point  between  the  membrane  and  the 
aluminum  mold.  A  third  solution,  considered  as  a  compromise,  was  to 
place  the  disc  in  an  epoxy  resin;  both  faces  of  the  membrane  were 
protected  by  a  piece  of  paper  tape,  prior  to  immersing  it  in  the 


. 


FIGURE  V-2 


IMBIBITION  CAPILLARY  PRESSURE  CELL 
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resin,  so  that  only  its  circumference  would  be  in  contact  with  the 
resin.  The  resin  imbibed  into  the  disc,  so  that  the  area  open  to 
flow  was  significantly  decreased.  It  was  finally  decided  to  glue  a 
previously  cleaned  membrane  in  a  titanium  ring  with  Devcon  F,  an  80% 
Aluminum  -  20%  Improved  Epoxy  Plastic.  Because  Devcon  F  cannot  be 
heated  to  high  temperatures,  the  glued  membrane  was  constantly  kept 
water-wet  by  letting  a  small  stream  of  distilled  water  flow  through  it. 
Thus,  the  wettability  was  not  allowed  to  change  significantly  between 
the  experiments.  This  method  was  found  to  be  fully  satisfactory. 

V. 3.2  Core  and  Cell  Preparation 

After  drainage  was  completed,  the  slices  to  be  used  in  the 
imbibition  experiment  were  rotated  at  maximum  speed,  in  order  to 
establish  a  saturation  as  close  to  irreducible  water  saturation  as 
possible.  The  slices,  kept  in  oil,  were  then  weighed.  The  cores 
were  reassembled  and  put  in  the  cells.  The  free  oil  and  water  levels 
were  very  carefully  measured  and  calculated  so  that  the  capillary 
pressure  at  the  outlet  end  of  the  core  would  be  close  to  zero,  but 
positive  in  order  to  prevent  the  water  from  flowing  out  of  the  core. 
Equilibrium  is  assumed  to  be  reached  when  the  free  water  level  does 
not  vary  for  a  sufficient  period  of  time. 

V. 3.3  Imbibition  Capillary  Pressure  and  Saturation  Measurement 

The  cell  was  designed  so  that  both  the  oil  and  water 
reservoirs  were  opened  to  the  atmosphere.  Under  this  condition,  the 
imbibition  capillary  pressure  was  expressed  as 
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 1.097  x  10 
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0) 


[D2(p  -p  )  +  r2p 
L  Mw'  vr w 


-  r2p  ] 
oMoJ 


(4-23) 


where 


r^  and  rQ  are  the  distances  from  the  rotary  axis  to  the 

free  water  and  oil  levels  respectively,  cm. 

D  is  the  distance  from  the  rotary  axis  to  the  point  where 

Pc  is  calculated,  cm. 

P  .  kPa. 
c 

The  average  saturation  in  each  slice  was  estimated  using  Eq . ( 5-2 ) : 


Equations  (4-23)  and  (5-2)  are  derived  in  Appendix  D. 


(5-2) 
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CHAPTER  VI 


RESULTS  AND  DISCUSSION 


VI. 1  DRAINAGE 

VI. 1.1  Experimental  Results 

The  samples  used  in  this  study  were  taken  from  three 
different  porous  rods  which  were  supposed  to  be  homogeneous.  Their 
physical  properties  are  given  in  Table  VI-1. 

The  fluids  used  were  always  water  and  kerosene;  this 
combination  was  chosen  because  lucite  is  quite  resistant  to  attack 
by  kerosene,  and  the  density  difference  between  both  fluids  is 
relatively  large,  which  is  favourable  to  attaining  high  capillary 
pressures.  The  properties  of  the  fluids  are  given  in  Table  VI-2. 

All  the  experiments  were  done  at  a  constant  temperature  of  about 
20°C,  and  the  fluid  properties  were  therefore  considered  to  be 
identical  for  all  the  runs. 

The  experiments  were  conducted  in  various  ways,  some  with 
slices  about  2.4  cm  long,  some  with  slices  about  1.2  cm  long.  Some 
runs  were  carried  out  at  one  single  speed  of  rotation;  in  other 
experiments,  the  speed  of  rotation  was  increased  by  steps,  and  the 
saturation  measured  after  each  step.  This  was  done  in  order  to  gain 
a  better  understanding  of  the  experimental  data  obtained,  because,  as 
it  has  already  been  pointed  out,  different  experimental  curves  can  be 
obtained  for  a  single  run  and  for  a  single  system  of  fluids  and 
porous  medium.  It  is,  therefore,  necessary  to  grasp  the  significance 


53 


54 


TABLE  VI-1 

CORE  PROPERTIES 


Sandstone 

P-6-C 

A1 undum 

P-12-C 

A1 undum 

Porosity,  % 

18.5 

36.9 

41.7 

Liquid  permeability,  ym2  0.131 

0.034 

0.155 

TABLE  VI-; 

2 

FLUID  PROPERTIES 

AT  20°C 

Densi ty 

Vi scosi ty 

kg/m3 

mPa-sec 

Water 

998.2 

1.00 

Kerosene 

788.3 

1.23 

Density  difference,  Ap 

=  209.9  kg/m3 

! 

Interfacial  tension,  y 

=  27.68  mN/m 

. 
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of  the  experimental  data,  and  to  be  able  to  distinguish  between  the 
capillary  properties  of  the  system,  and  the  effects  due  to  the 
parameters  of  the  experiment,  such  as  length  of  the  slices  and  speed 
of  rotation.  These  parameters  are  tabulated  in  Table  E-l  through 
E-8  in  Appendix  E. 

No  problem  was  encountered  with  the  centrifuge,  whose 
flexible  drive  shaft  and  feedback  speed  control  allowed  the  rotation 
to  be  smooth  and  steady. 

Experimental  curves  are  shown  in  Figure  VI - 1  to  Figure  VI-8, 
and  data  are  given  in  Table  F-l  through  F-8  in  Appendix  F. 


VI. 1.2  Error  Analysis 

The  drainage  capillary  pressure  is  given  by  Eq. (4-9) . 
Therefore,  the  relative  error  on  capillary  pressure  is  given  by 


_  .  ,  6h 

{PC  «(Ap)  .  ,  x  Sre  2  A  6h 

T7  =  ~  r  h  h 

re  2 

where 

6  (ap)  ^  5  x  10'4  g/cm3 

6a )  -  10  rpm  at  1000  rpm 

-2 

6r  -  10  cm 

e 

6  h  —  10"3  cm 


Thus,  the  relative  error  for  capillary  pressure  should  not  exceed  3 I. 
The  average  saturation  was  measured  by  Eq. (5-2) .  The 


relative  error 


is, 


therefore, 

6M 1  +  6M, 
w  l 
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FIGURE  V 
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where 

6M 1  ^  5  x  10"4  g 

w  v 

6M.  ^  5  x  10“4  g 

«Vp  =S  10'3  cm3 

<5(Ap)  =:  5  x  10'4  g/cm3 

Hence,  the  relative  error  for  the  measured  average  saturation  should 
not  be  more  than  about  3%. 

VI. 1.3  Numerical  Results 

The  Marquardt  algorithm  for  least-squares  estimation  of 

nonlinear  parameters  (96)  was  used  to  fit  the  data  points  by 

estimating  the  parameters  which  define  the  model.  Capillary  pressure 

was  the  independent  variable,  and  S~  was  the  dependent  variable, 

w 

according  to  Eq . (4-17)  and  Eq . (4-18). 

The  calculated  values  of  S~  are  given  in  Table  F-l  to 

w 

Table  F-8  of  Appendix  F,  and  the  generated  capillary  pressure- 
saturation  curves  are  shown  in  Figure  VI-9  to  Figure  VI-11.  The 
estimated  parameters  are  given  in  Table  VI-3,  along  with  the 
associated  linear  confidence  limits. 

VI. 2  IMBIBITION 

Imbibition  experiments  were  conducted  on  cores  made  of  four 
slices;  each  slice  was  about  1.2  cm  long.  It  was  not  possible  to  use 
six  slices,  as  in  the  drainage  experiments,  because  it  was  necessary 
to  keep  the  capillary  pressure  at  the  bottom  of  the  cores  positive. 
Prior  to  starting  the  experiments,  the  irreducible  water  saturation 


;  ■ 


Capillary  Pressure, 


66 


Capillary  Pressure, 


67 


TABLE  VI-3 


68 


OO 


OO 

LU 

Cd 


<c 

O 


i — i 


Cd 

LU 


L- 

<U  0J 


"O  > 

CO 

LO 

CO 

cxj 

CO 

o 

o 

cr  s- 

CO 

o 

r— 

cr. 

r— 

CO 

co 

ID  D 

CO 

CO 

o 

00 

co 

CO 

LO 

LO 

CD 

CO 

r^. 

o 

co 

r — > 

CXI 

co 

CXI 

, - 

<U  CU 

1— 

CO 

CXJ 

CO 

CO 

CO 

LO 

LO 

S_  -d 

<  +-> 

o 

C\J 

CXI 

CXI 

CXI 

CXI 

CXI 

CXJ 

CXJ 

S-  CU 

i 

1 

1 

| 

1 

1 

1 

S-  4-> 

o 

o 

o 

o 

o 

o 

o 

o 

LU  fO 

r— 

r— 

r— 

r— 

r— 

f— 

r— - 

f— 

E 

T3  -i- 

X 

X 

X 

X 

X 

X 

X 

X 

S-  -U> 

fO  CO 

o 

CO 

LO 

cr. 

CXI 

LO 

CXJ 

“D  LU 

CO 

cr. 

CVJ 

CO 

CO 

CO 

c 

o 

r— 

o 

co 

oo 

CXJ 

fO  4- 

4->  O 

T— 

i — 

r— 

o* 

o 

o 

1— 

1— 

00 

o 

o 

o 

CXJ 

CO 

LO 

c 

o 

o 

o 

cr. 

r^ 

r>. 

CO 

o 

c 

1 

1 

1 

o 

o 

CD 

o 

1 

LO 

co 

co 

o 

00 

«vT 

f— 

r— 

cr. 

r— 

CO 

LO 

CO 

LO 

CO 

i^> 

cr. 

CO 

o’ 

o 

o 

o 

CD 

o 

r— 

*o 

+1 

+1 

+! 

+1 

+i 

+1 

+i 

+1 

D- 

r-- 

r— 

CO 

r-'. 

r— 

CXI 

CTi 

CXI 

LO 

o 

cr. 

cr, 

cxj 

CXJ 

CXJ 

CO 

CO 

r-^ 

cr. 

LO 

CO 

• 

• 

• 

• 

• 

• 

• 

• 

CXI 

cr i 

o 

co 

CXJ 

CXJ 

CXJ 

CO 

1 

i 

CXI 

CO 

o 

LO 

CO 

■^r 

cn 

CXI 

CXI 

co 

r— 

r— 

r— 

CXJ 

r-* 

CO 

LO 

1 — 

o 

CXI 

• 

• 

• 

• 

• 

• 

• 

• 

CXI 

1— 

*— 

•— 

CXJ 

•— 

CO 

fc> 

+1 

+1 

+1 

+i 

+1 

+i 

+1 

+i 

CD 

o 

, _ 

o 

CO 

oo 

CO 

CXI 

LO 

CO 

<XJ 

r^ 

o 

o 

CO 

CO 

00 

1 - 

LO 

CXI 

• 

• 

• 

• 

• 

• 

• 

• 

i_n 

o 

00 

CO 

o 

*=|- 

LO 

1 

CXI 

1 

CO 

CO 

CO 

CO 

CXI 

cn 

CXJ 

cn 

CO 

CO 

, — 

o 

r— 

o 

CXI 

CXJ 

o 

o 

o 

o 

c 

o 

o 

o 

• 

• 

• 

• 

• 

• 

• 

o 

o 

o 

o 

o 

o 

o 

o 

•r— 

+1 

+1 

+1 

+1 

+1 

+1 

+1 

+1 

00 

LD 

CO 

CO 

, _ 

CO 

cr. 

CXI 

CO 

r^. 

CXI 

LO 

r^. 

CO 

1 - 

o 

o 

o 

1 - 

' — 

• 

• 

• 

• 

• 

• 

• 

o 

o 

o 

o 

o 

o 

o 

o 

c 

Z3 

CXI 

CO 

LO 

CO 

r-^ 

CO 

QC 

69 


was  established  in  each  slice.  When  the  saturation  of  the  inner  slice 
was  S  • ,  this  slice  was  removed,  and  the  core  was  raised.  The  second 

W  I 

slice  of  the  original  core  was  then  the  top  slice,  and  the  bottom 

slice  was  replaced  by  an  end  butt  made  of  the  same  material,  and  100% 

water  saturated.  The  centrifuge  was  restarted,  and  S  .  was 

wi 

established  in  the  top  slice;  this  procedure  was  repeated  until  S  . 
was  obtained  in  eight  slices,  that  is  in  two  complete  cores. 

Extreme  care  was  taken  in  the  measurement  of  the  water  and 
oil  levels,  in  their  respective  reservoirs.  The  duration  of  each 
experiment  was  about  one  week;  the  last  experiment,  conducted  on 
P-12-C,  was  carried  on  during  24  days,  but  in  spite  of  the  fact  that 
the  capillary  pressure  at  the  inner  end  of  the  core  was  kept  as  low  as 
6.79  kPa,  no  imbibition  whatsoever  was  obtained. 

VI. 3  DISCUSSION 

VI. 3.1  Drainage 

VI. 3. 1.1  Experimental  Results 

Run  2  and  Run  7  were  conducted  on  cores  made  of  three  slices 
about  2.4  cm  long;  three  speeds  of  rotation,  in  the  range  1000  to 
3000  rpm,  were  chosen  in  both  runs,  in  order  to  obtain  the  complete 
capillary  pressure  curves.  The  saturation,  S^,  in  each  slice,  was 
calculated  after  equilibrium  was  obtained  at  a  particular  speed,  as 
well  as  the  capillary  pressure,  P^,  at  the  top  of  each  slice,  for 
that  particular  speed.  The  raw  data  thus  obtained  for  Run  7  are 
plotted  in  Figure  VI-12.  Data  obtained  for  Run  2  would  lead  to  a 
similar  kind  of  plot.  Another  interpretation  of  the  data  must 
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therefore  be  resorted  to,  in  order  to  get  a  more  conventional  type 

of  experimental  curve.  Szabo  (87)  obtained  his  curves  by  plotting 

the  calculated  capillary  pressure,  P  ,  at  the  center  of  each  slice, 

against  the  related  average  saturation,  S~.  Figure  VI-13  shows  the 

w 

curve  thus  obtained  for  Run  7.  The  improvement  is  evident,  but  the 
saturation  in  the  lower  part  of  the  curve  is  still  double-valued. 


Such  behaviour  may  be  easily  understood.  Figure  VI-14 


shows  a  slice  and  the  saturation  profile  within  it.  It  is  evident 

that  a  plot  of  P  versus  S~  is  nearer  to  a  true  capillary  pressure- 

c  c  w 

saturation  curve,  than  is  a  plot  of  P  .  against  S~.  However,  there 

Cl  w 

is,  a  priori,  no  reason  why  Pcc  should  correspond  exactly  to  Sw,  and 

a  plot  of  P  against  ~  cannot  be  accepted  as  a  true  capillary 
r  cc  w 

pressure-saturation  curve. 

Figure  VI-15  shows  two  similar  cores,  one  spinning  at  mi 
rpm  and  the  other  spinning  at  m2  rpm.  The  values  of  Pcl-  at  the  top 
of  each  slice  are  also  indicated.  Let  m2  be  such  that  p2  =  p3- 
By  vi rtue  of  Eq. (4-9) , 


k  =  3h 


(6-1) 


and 


£  =  2h 


(6-2) 


Therefore 


and 


(6-3) 


►  f 
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FIGURE  VI-14  GRAVITY  EFFECT 


FIGURE  VI-15  EFFECT  OF  SPEED  VARIATION,  u>!  <  w2 
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It  follows  that 

Pi  =  AApc4  (rg  -  j)  h 

2  3(2r  -  3h) 

=  AAp“i  '  4(r  -  h)  (re  '  P  h  (6-4) 

e 

A1  so, 

P2  =  AApwi  (re  -  h)  2h  (6-5) 

By  comparing  Eq . (6-4)  and  Eq.(6-5),  it  is  easy  to  see  that 

P2  >  Pi  (6-6) 

Consequently 

S3  <  S2  (6-7) 

Therefore,  it  is  not  possible  to  obtain  a  single  curve  by  plotting 
the  raw  data,  when  different  speeds  of  rotation  are  used  in  the 
course  of  one  experiment. 

Another  interpretation  of  the  data  was  necessary;  that  is, 
the  average  saturation,  ST",  was  calculated  from  the  outlet  end  of 
the  core,  to  the  plane  corresponding  to  P  . .  Such  curves  were 

I 

plotted  for  Run  2  and  Run  7  [Figure  VI-2  and  Figure  VI-7].  They  are 
much  more  satisfactory. 

Let  Pi  =  P2;  then,  according  to  Eq . (4-1 2 )  and  Eq . (4-1 7) , 
and  combining  Eq . (4-1 2 )  with  Eq.(6-3), 
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Hence,  in  spite  of  the  fact  that  =  P2,  sf  is  not  strictly  equal 
to  S^.  Figures  VI-16  through  VI-18  show  that  the  corrected  data 
pertaining  to  the  bottom  slices  are  not  on  the  modelled  curve. 


This  system  is,  therefore,  peculiar  in  that  it  does  not 
allow  the  obtainment  of  conventional  experimental  curves,  if 
different  speeds  of  rotation  are  used  in  a  single  experiment. 
Moreover,  the  time  required  to  obtain  equilibrium  at  each  speed  of 
rotation  was  found  to  be  quite  long,  up  to  47  hours  for  the  tightest 
core. 

The  cores  used  in  Run  1,  Run  5,  Run  6,  and  Run  8  were  made 
of  six  slices,  each  about  1.2  cm  long.  Increasing  the  number  of 
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slices  has  the  advantage  that  more  points  are  obtained  at  once,  and, 
if  the  size  of  the  slices  and  speed  of  rotation  are  adequately  chosen, 
an  entire  capillary  pressure-saturation  curve  can  be  obtained  at  one 
speed  of  rotation.  Because  the  bottom  slice  was  always  partially 
desaturated  after  each  run,  the  experimental  capillary  displacement 
pressure  was  obtained  by  extrapolation;  this  measurement  cannot, 
therefore,  be  considered  as  accurate.  On  the  other  hand,  the 
irreducible  water  saturation  was  very  nearly  obtained.  Again,  two 
experimental  curves  can  be  drawn  by  plotting  either  P  .  against  Sw 
or  P  .  against  IT"  [Figure  VI-19  and  VI-6].  Neither  of  these  curves 
represents  the  true  capillary  pressure  curve. 

In  order  to  get  a  better  view  of  what  the  true  capillary 
pressure  curve  should  look  like,  Run  3  and  Run  4  were  conducted. 

They  are  similar  to  Run  2  and  Run  5  respectively,  with,  however,  the 
following  modification:  the  bottom  of  each  core  was  cut  in  two  slices; 
the  outer  slice  was  about  0.5  cm  long,  and  the  effect  of  gravity  is, 
therefore,  expected  to  become  small  due  to  a  small  variation  of  Pc 
from  the  inlet  to  the  outlet  face  of  the  slice.  It  is  very  interesting 
to  note  that: 

i)  the  points  corresponding  to  the  thin  outer  slices  do 
not  fall  on  the  experimental  capillary  pressure  curve  [Figure  VI-3], 

ii)  although  they  display  an  inflection  point  near 
maximum  saturation,  these  data  are  much  closer  to  the  model  than  the 
data  obtained  for  the  large  slices  {Figure  VI-20], 

iii)  the  extrapolated  capillary  displacement  pressure  is 
drastically  different  and  smaller  than  that  obtained  by  extrapolating 
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the  data  pertaining  to  the  large  slices  [Figure  VI-3  and  Figure  VI-20]. 

iv)  the  saturation  measured  in  these  small  slices  is  still 
an  average  saturation  and  is,  therefore,  slightly  too  high.  A  true 
capillary  pressure  curve  should,  therefore,  slightly  deviate  from 
these  points  and  be  located  to  their  left  [Figure  VI-20]. 

VI. 3.1 .2  Numerical  Results 

The  program  used  to  estimate  the  parameters  makes  it 
possible  to  let  all  the  parameters  float,  or,  conversely,  to  fix  any 
chosen  parameter.  If  some  of  these  parameters  are  strongly  related 
to  each  other,  a  variation  of  one  of  them  induces  a  variation  in  the 
others,  and  conversely.  In  such  a  case  the  algorithm  diverges  very 
quickly.  This  was  the  case  for  n  and  a,  which  determine  the  area 
under  the  curve  and  are  very  strongly  related;  this  area  is  given 
by 

area  =  +  Pd  (4-7) 

The  term  is  represented  in  Figure  VI-21  by  the  shaded  portion 
of  the  area  under  the  curve.  Therefore,  it  was  necessary  to  fix  n. 

The  three  other  parameters,  S  . ,  a,  and  P^,  were  allowed  to  float, 
and  convergence  was  obtained  in  all  cases  after  five  to  nine 

_5 

iterations,  with  a  convergence  criterion  e  =  10  .  Another  value 

of  n  was  then  chosen,  and  the  program  rerun.  The  standard  error  of 
estimate  was  the  criterion  used  to  determine  which  value  of  n  gave 
the  best  fit.  In  all  but  one  run,  Run  4,  the  standard  error  of 
estimate  was  found  to  increase  continuously  with  n,  or,  on  the 
contrary,  to  decrease  continuously  with  an  increasing  n.  In  the  first 
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case,  n  was  lowered  until  a  mathematical  impossibility  would  obtain; 
the  last  value  of  n  for  which  the  calculations  could  be  carried 
through,  and  the  corresponding  values  of  S  .,  P,  and  a,  were  chosen 
as  the  best  estimates  of  the  parameters  for  the  particular  run.  In 
such  cases,  the  logarithmic  model  was  found  at  times  to  be  totally 
inadequate,  because  no  combination  of  initial  parameters  could  be 
found  that  would  allow  the  program  to  converge.  On  the  other  hand, 
when  the  standard  error  of  estimate  was  found  to  decrease  for 
increasing  n,  the  logarithmic  model  converged  easily  and  gave  the 
best  fit;  it  was  not  possible  to  increase  n  continuously  to  one  due 
to  prohibitively  high  values  of  the  exponent  y~,  involved  in  the 
calculations,  when  n  approaches  one.  The  highest  acceptable  value 
of  n  was  found  to  be  of  the  order  of  0.97,  after  which  the  logarithmic 
model  had  to  be  used. 

_2 

The  standard  error  of  estimate  ranged  between  0.484  x  10 

p  _2 

and  1.835  x  10"  ,  which  is  quite  similar  to  the  values  0.368  x  10 
and  1.537  x  10"2,  obtained  by  An! i  (97).  The  linear  confidence 

limits  for  the  parameters  are  given  in  Table  VI-3;  they  are 
consistently  higher  for  a  than  for  S^-  and  P^,  a  fact  which  has  been 

already  pointed  out  by  An 1 i  (98). 

It  seems  that  the  estimates  of  the  parameters  are  fairly 
constant  for  comparable  experiments,  but  that  they  vary  widely  with 
the  experimental  procedure,  that  is  the  size  and  number  of  slices  and 
speed(s)  of  rotation.  Hence,  they  are  comparable  for  Run  2  and  Run  3 
on  the  one  hand,  and  comparable  for  Run  5  and  Run  6  on  the  other 
hand,  but  the  two  sets  thus  obtained  are  very  different.  This  is  not 


surprising,  however,  because  the  model  is  not  grounded  on  physical 
considerations  but  is  merely  a  fit;  therefore,  it  is  not  reasonable 
to  expect  the  individual  parameters  to  remain  constant,  even  though 
the  conditions  of  the  experiments  may  change.  This  is  true 
particularly  in  the  case  where  the  parameters  are  obtained  by 
extrapolation.  Run  4  is  an  illustration  of  this  fact:  the  estimated 
value  of  Pd  for  this  run  is  negative,  in  spite  of  the  fact  that  the 
core  was  water-wet. 

Strictly  speaking,  a  model  which  does  not  exactly  reflect 
the  physics  of  the  phenomenon  under  study  should  not  be  extrapolated 
beyond  the  range  of  available  experimental  data,  because  every  point 
has  an  equal  weight  in  the  determination  of  the  parameters;  the  best 
estimates  of  the  latter  are,  therefore,  best  only  for  the  given  set 
of  experimental  data.  However,  if  a  combination  of  two  or  more 
parameters  can  be  given  a  physical  meaning,  it  is  worthwhile 
investigating  this  combination. 

It  has  been  previously  stated  that  the  area  under  the 
capillary  pressure-saturation  curve  corresponds  to  the  reversible 
external  work,  and  if  equilibrium  has  been  reached,  one  may  expect  th 
work  to  be  almost  constant  and  independent  of  the  experimental 
procedure.  The  value  of  the  area  under  each  curve  is  given  in 
Table  VI-3.  The  average  for  sample  P-12-C  is  31.29  with  a  standard 
deviation  of  2.20;  for  P-6-C,  the  average  is  52.06  with  a  standard 
deviation  of  0.75.  These  values  are,  therefore,  meaningful.  Of 
particular  interest  is  Run  4;  the  value  of  the  area  under  this  curve 
falls  within  2 %  of  the  average  obtained  for  sample  P-12-C,  although 
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the  estimated  parameters,  for  this  run,  are  particularly  far  from 
those  obtained  in  the  other  runs,  and  even  meaningless  as  far  as 
is  concerned. 

VI. 3. 1.3  Capillary  Displacement  Pressure  -  End  Effect 

It  is  commonly  accepted  that  extrapolating  experimental 
data  to  Sw  =  1  yields  P^.  Anli  (99)  reports  the  observation  of 
a  significant  discrepancy  between  the  values  of  estimated  and 
extrapolated  capillary  displacement  pressure  for  all  his  runs  and 
suggests  that  there  is  a  linear  correlation  between  the  discrepancy 
and  the  permeability  of  the  cores;  the  discrepancy  increases  with 
permeability.  The  largest  discrepancy  was  obtained  for  his  shortest 
core  and  does  not  fit  the  correlation;  neither  does  the  discrepancy 
for  the  tightest  core. 

Anli  (100)  points  out  that  the  displacement  capillary  pressure, 
obtained  by  extrapolation,  is  expected  to  be  higher  than  the  true 
capillary  displacement  pressure  and  argues  that 

"the  deviation  of  the  estimated  Pd  from  that  observed 
is  not  necessarily  a  quantitative  representation  of 
the  end  effect,  but  could  also  be  a  result  of  the 
heterogeneity  of  the  sample". 

It  would,  indeed,  be  difficult  to  estimate  quantitatively 
the  end  effect  or  the  role  played  by  gravity  and  sample  heterogenei ty. 
But  one  can  argue  that  if  the  discrepancy  between  estimated  and 
extrapolated  were  due  to  heterogeneity  of  the  samples,  there  would 
be  no  apparent  reason  for  which  this  discrepancy  always  goes  in  the 
same  direction.  According  to  Szabo  (89),  the  measured  capillary 
pressure  curve  is  always  under  the  true  capillary  pressure 
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curve;  therefore,  the  discrepancy  would  always  be  opposite  to  that 
obtained  by  An! i  (99).  The  reason  is  that  Szabo  (89)  plots  the 
capillary  pressure  at  the  bottom  of  the  core  against  S  ,  whereas 
Anli  (101)  plots  the  capillary  pressure  at  the  top  of  the  core  against 
Sw.  The  true  capillary  pressure  curve  is,  of  course,  expected  to  be 
somewhere  in  between.  The  deviation  is,  therefore,  more  likely  due 
to  a  combination  of  gravity  effect  and  end  effect,  than  to  the 
heterogeneity  of  the  samples. 

It  is  worthwhile,  at  this  point,  to  investigate  Eq . ( 4-9 ) 
in  connection  with  the  end  effect.  In  order  for  the  water  to  be 
expelled  from  a  water-wet  porous  medium,  there  must  be  a  capillary 
pressure  inversion  inside  the  core,  so  that  the  pressure  in  the  water 
phase,  at  the  outlet  end  of  the  core,  is  higher  than  the  pressure  in 
the  oil  phase  at  the  same  point  (102).  Therefore,  the  assumption 
commonly  accepted,  that  P  =  0  at  the  outlet  end  of  the  core,  is 
always  violated.  This  results  in  higher  displacement  and  capillary 
pressures  than  what  would  be  obtained  without  end  effect,  until  the 
end  effect  becomes  so  small  that  it  can  be  neglected. 

The  drainage  capillary  pressure  equation  is 

P  =  AApoo2  (r  -  j)  h  (4-9) 

G  Q  £ 

and  the  displacement  capillary  pressure  is  given  by 

Pd  =  AaPco2  (re  -  *£•)  h* 


where  h*  is  the  location  where  Pc  =  P^  for  the  angular  speed  to. 
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Pd  has  a  finite,  positive  value  for  a  water-wet  core. 

Therefore,  when  to  goes  to  infinity,  h*  must  go  to  zero,  and  Pd  is 

finite,  not  equal  to  zero.  But,  by  hypothesis,  P  =  0  at  h  =  0. 

0 

Thus  there  is  a  contradiction. 

Let  h°(oo)  be  the  location  at  which  P  =  0.  One  has: 

c 

h  +  h° 

Pc  =  Aapu2  (re - ^ - )<h  -  h°)  (6-10) 

Then,  at  the  outlet  end  of  the  core, 

h° 

Pc(h  =  0)  =  -  AApu2  (re  -  -j-)  h°  (6-11) 

When  w  goes  to  infinity,  and  for  P  (h  =  0)  to  be  finite,  it  is 
necessary  that 

1  im  h°  =  h0(°°)  =  0 

This  was  expected  a  priori,  because  the  end  effect  becomes  very  small 
at  high  speeds  of  rotation.  Therefore,  if  an  experiment  is  conducted 
at  a  single  speed  of  rotation,  as  was  the  case  in  Runs  1,  4,  5,  6  and 
8,  h°  is  fixed,  whereas  it  varies  in  experiments  such  as  Runs  2,  3 
and  7.  In  the  first  case,  if  P  and  PQ2  are  the  capillary  pressures, 
calculated  with  Eq . (6-10),  at  hi  and  h2, 

h  i  +  h2 

Pci  -  PC2  =  AaP0)2  (re - 2 - )(hl  '  h2}  (6“12) 

The  difference  in  capillary  pressures,  calculated  at  two  different 
locations  within  the  core,  is  independent  of  h  .  Therefore,  the  end 
effect  is  only  to  shift  the  curve  upwards.  Conversely,  the  capillary 
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pressure  at  the  same  location  in  a  core  rotated  at  two  different 
speeds,  wj  and  u)2,  is  dependent  upon  h?  and  h2,  where  h?  and  h2  are 
the  locations  at  which  P  =  0  for  the  angular  speeds  and  w2 

V 

respectively: 

h°  h° 

PC1-  PC2  ■  AAp[(re  -  |)h(^  -  4)  -  u*h?(re  -  -2.)  +  “2h°(^e  -  T)] 

(6-13) 

Because  P  increases  with  to,  whereas  h°  decreases  with  to,  the 
discrepancy  due  to  the  end  effect  is  expected  to  increase  towards 
high  saturations  and  to  be  the  highest  at  S  =  1;  that  is,  for  the 
displacement  capillary  pressure.  Therefore,  an  experimental 
procedure  whereby  the  angular  speed  is  varied  to  obtain  the 
capillary  pressure  curve  is  more  subject  to  errors  due  to  the  end 
effect,  than  is  a  method  whereby  the  curve  can  be  obtained  at  a 
single  high  speed  of  rotation. 

VI. 3. 2  Imbibition 

A  rigorous  analysis  of  the  system  under  study  is  difficult. 
Due  to  the  difference  in  pressure  gradients  in  the  oil  and  water 
phases,  the  capillary  pressure  was  always  higher  at  the  inner  end  of 
the  core  than  at  the  outer  end.  The  semi -permeable  membrane,  placed 
between  the  water  reservoir  and  the  inner  end  of  the  core,  was  fully 
water  saturated,  which  resulted  in  a  saturation  discontinuity  at  the 
contact  plane  between  the  core  and  the  membrane.  The  core  and 
membrane  were  kept  in  capillary  contact  by  means  of  a  disc  of 
Kleenex  paper. 
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If,  for  a  given  r  and  r  ,  the  speed  of  rotation  is  chosen 

w  u 

so  that  the  capillary  pressure  at  the  inner  end  of  the  core, 

calculated  with  Eq . ( 4-23 ) ,  corresponds  to  S  . ,  imbibition  is  not 

w  i 

possible,  and  no  change  occurs  in  the  system.  If  the  speed  of 
rotation  is  decreased  sufficiently,  so  that  the  forces  due  to 
interfacial  tensions  and  wettability  overcome  the  force  due  to  the 
excess  pressure  in  the  oil  phase  as  compared  to  the  pressure  in  the 
water  phase,  imbibition  should  occur. 

At  the  start  of  the  experiment,  the  capillary  pressure 
inside  the  core  is  much  larger  than  the  capillary  pressure  calculated 
with  Eq . ( 4-23 ) .  If  imbibition  occurs,  rw  increases  and  rQ  decreases; 

therefore,  the  calculated  capillary  pressure  increases  in  the  course 
of  the  experiment,  whereas  it  decreases  in  the  core  as  a  result  of 
imbibition.  It  is  therefore  believed  that  the  water  should  progress 
towards  the  outer  end  of  the  core,  and  that  the  location  at  which 
equilibrium  obtains  should,  in  the  same  manner,  progress  from  the 
inlet  to  the  outlet.  Let  us,  indeed,  suppose  that  at  a  time  Tls  a 
state  of  equilibrium  is  obtained  at  the  distance  from  the  inlet 
end  of  the  core;  if  the  state  of  equilibrium  progresses  outwardly, 
at  time  Tj.  +  <ST,  the  imposed  hydrostatic  pressure  has  increased  at 
the  location  £1}  as  a  result  of  imbibition  in  the  remainder  of  the 
core.  Consequently,  the  saturation  at  £}  must  decrease  by  a 
corresponding  amount,  in  order  for  a  new  state  of  equilibrium  to 
obtain.  This  is  possible  because,  by  hypothesis,  the  remainder  of 
the  core  is  not  yet  at  equilibrium,  and  the  excess  water  at  £i  can 
move  towards  the  outer  end  of  the  core.  On  the  contrary,  if 
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equilibrium  were  first  obtained  at  the  outer  end  and  progressed 
towards  the  inner  end  of  the  core,  there  would  be  no  way  for  the 
excess  water  in  the  lower  part  of  the  core  to  escape,  as  the  imposed 
hydrostatic  pressure  increases  as  a  result  of  further  imbibition  in 
the  upper  part  of  the  core. 

As  a  result  of  this  analysis,  it  seems  that  the  saturation 
at  one  point  inside  the  core  oscillates  continuously,  until 
equilibrium  is  eventually  reached  in  the  entire  core.  This  is  not 
desirable,  because  it  is  equivalent  to  an  alternation  of  imbibition 
and  drainage  cycles.  The  amplitude  of  the  oscillations  is  not 
measurable,  and,  as  a  result,  it  is  not  possible  to  evaluate  the 
consequences  of  this  effect  on  the  experimental  results. 

Another  drawback  of  the  method  is  that,  whereas  it  was 
desirable  to  have  a  high  Ap  in  the  drainage  experiments,  it  is 
desirable  to  have  a  small  Ap  in  the  imbibition  experiments.  Indeed, 
if  L  is  the  length  of  the  core,  and  if  PQ  and  Pci  are  the  capillary 
pressures  at  the  outer  and  inner  ends  respectively  of  the  core,  then 

Pci  =  Po  *  Tl  (2re  -  L)  Ap  (6-14) 

In  order  for  imbibition  to  occur  quickly,  it  is  desirable 

to  use  as  high  an  w  as  possible,  and  to  keep,  at  the  same  time,  a 

fairlv  low  P  •;  therefore  Ap  should  be  small.  This  would  be  even 
J  ci 

more  critical  in  the  case  of  unconsolidated  porous  media,  for  which 
the  capillary  pressures  involved  are  small. 

The  reasons  for  the  failure  of  these  experiments  are  not 

known  with  certainty.  The  cores  were  water-wet.  Therefore,  the 


most  probable  reason  is  that  the  capillary  pressure  at  the  inner  end 
of  the  core  was  too  high.  This  suggests  that  it  should  be  extremely 
small;  the  minimum  possible  speed  of  rotation  of  the  centrifuge  was 
500  rpm;  therefore,  a  lowering  of  capillary  pressure  would  be  possibl 
only  by  controlling  r  and  r  .  But  this  control  is  not  sufficient  at 

W  0 

very  low  capillary  pressures.  Indeed,  when  the  term 
[(re  -  h)2(po  -  pw)  +  r2  pw  -  r2  pQ]  of  Eq . (4-23)  becomes  very  small, 
a  small  error  in  the  measurement  of  r  has  an  important  effect,  and 
the  term  may  easily  become  negative.  Moreover,  the  centrifugal 
acceleration  distorts  somewhat  the  free  oil  and  water  levels. 
Therefore,  the  calculated  values  of  P  would  be  totally  erratic.  In 
his  experiments,  Szabo  (103)  was  able  to  use  speeds  as  low  as  195 
rpm  and  had,  therefore,  better  control  over  the  value  of  P  . 


CHAPTER  VII 


CONCLUSIONS  AND  RECOMMENDATIONS 

An  attempt  has  been  undertaken  to  gain  a  better 
understanding  of  the  significance  of  the  experimental  capillary 
pressure  -  average  saturation  data  obtained  with  the  centrifuge 
method  on  multi-slice  cores.  The  size  and  number  of  slices  were 
varied  in  the  experiments,  and  one  or  several  speeds  of  rotation  were 
used  in  the  course  of  the  experiments.  A  new  mathematical  model, 
to  fit  the  data,  was  studied  and  compared  with  another  model  already 
studied  in  the  literature.  Attempts  have  been  undertaken  to  obtain 
the  imbibition  branch  of  the  hysteresis  loop  with  the  centrifuge. 

As  a  result  of  this  investigation,  the  following  conclusions  can  be 
drawn;  it  is  pointed  out  that  conclusions  1,  2  and  4  are  relevant 
to  the  centrifuge  method  conducted  on  multi-slice  cores,  whereas 
conclusions  3,  5,  6  and  7  are  relevant  to  both  the  conventional 
centrifuge  method  and  the  quick  method: 

1)  the  experimental  data  vary  with  the  parameters  of  the 
experiments,  such  as  size  of  the  slices  and  whether  one  or  several 
speeds  of  rotation  are  used  in  the  course  of  the  experiments. 

2)  at  least  two  experimental  curves  can  be  obtained  for 
a  given  set  of  data,  depending  on  whether  capillary  pressure  is 
plotted  against  the  average  saturation  in  each  slice  or  against  the 
average  saturation  between  the  outer  end  of  the  core  and  the  plane 

at  which  P  is  calculated.  This  suggests  that  experimental  data 
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alone  do  not  represent  the  true  capillary  properties  of  the  medium 
under  study. 

3)  extrapolated  data,  such  as  displacement  capillary 
pressure,  are  not  reliable,  especially  when  the  size  of  the  slices  is 
large  enough  for  the  gravity  effect  to  be  important.  This  was  already 
shown  by  Anli. 

4)  if  different  speeds  of  rotation  are  used  in  the  course 
of  an  experiment  conducted  on  a  multi-slice  core,  a  conventional 
capillary  pressure  -  average  saturation  curve  cannot  be  obtained 
experimental ly. 

5)  the  assumption  that  P  =  0  at  the  bottom  of  the  core, 
although  it  is  probably  acceptable  in  view  of  the  accuracy  of  the 
data  obtained,  bears  in  itself  a  contradiction;  this  contradiction 
can  be  removed  only  by  introducing  the  notion  of  end  effect. 

6)  although  the  end  effect  is  probably  small,  data  obtained 
by  varying  the  speed  of  rotation  are  more  likely  to  be  distorted  than 

t 

those  obtained  at  a  single  high  speed  of  rotation. 

7)  the  most  probable  reason  for  the  failure  of  the 
imbibition  experiments  is  that  the  capillary  pressure  at  the  top  of 
the  core  could  not  be  sufficiently  lowered.  The  centrifuge  method, 
however,  does  not  seem  to  be  a  very  sound  method  for  obtaining 
imbibition  data,  because  water  must  penetrate  the  core  on  the  side 
where  the  saturation  is  the  lowest;  that  is,  on  the  side  where  the 
permeability  to  water  is  the  lowest. 
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The  following  suggestions  can  be  made: 

1)  drainage  data  could  be  improved  by  using  slightly 
larger  cores,  made  of  slices  about  0.7  cm  long,  so  that  the  outer 
slice  would  remain  almost  fully  saturated,  whereas  the  inner  slice 
would  have  irreducible  water  saturation,  when  using  a  single,  well 
chosen,  speed  of  rotation. 

2)  displacement  capillary  pressure  should  be  measured  by 
the  diaphragm  method,  because  this  method  is  free  from  end  effect. 

3)  once  and  S  .  are  known  experimental ly,  the  best 
values  of  n  and  a  could  be  evaluated  by  using  the  proposed  model. 

The  curve  thus  generated  is  expected  to  give  a  fairly  good 
representation  of  the  true  capillary  properties  of  the  system  under 
study. 

4)  if  it  is  desired  to  have  a  single  value  of  n  for  the 
drainage  model,  a  value  of  0.90  would  be  suggested;  this  is  not  based 
on  theoretical  considerations,  but  it  was  found  that  with  such  a 
value  of  n,  it  was  always  easy  to  find  initial  values  of  the 
parameters  which  allowed  the  algorithm  to  converge.  Moreover,  the 
variation  of  the  standard  error  of  estimate  with  n  is  very  small,  and 
the  values  of  the  parameters  thus  estimated  are  close  to  those  which 
define  the  best  estimate. 

5)  unless  a  centrifuge  capable  of  very  low  speeds  of 
rotation  is  available,  it  does  not  seem  advisable  to  use  the 
centrifuge  method  on  unconsolidated  cores,  or  to  try  to  obtain 
imbibition  curves.  If  imbibition  data  are  obtained  with  the 
centrifuge,  it  would  be  necessary  to  compare  the  results  to  those 
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obtained  with  the  diaphragm  method  on  a  similar  system  of  fluids  and 
porous  medium,  in  order  to  check  the  validity  of  the  centrifuge  method. 

6)  a  method  could  possibly  be  devised,  so  as  to  yield  the 
true  capillary  pressure  curve  of  a  system  by  using  dye  and  measuring 
the  density  of  light  passing  through  transparent  pictures  of  the  faces 
of  the  slices.  Because  the  amount  of  dye  is  proportional  to  the 
amount  of  water,  there  is  a  relationship  between  the  amount  of  light 
transmitted  and  the  saturation. 

7)  care  should  be  taken  in  the  core  preparation.  If  the 
cores  contain  large  amounts  of  water  of  crystallization,  they  should 
not  be  heated  up  to  high  temperatures  in  order  to  prevent  evaporation 
and  subsequent  uptake  of  water  of  crystallization  during  the 
saturation  process. 

8)  although  no  problem  of  false  minima  was  encountered  in 
the  course  of  this  study,  this  has  not  always  been  the  case. 
Consequently,  great  care  should  be  taken  in  using  the  Marquardt 
algorithm.  It  is  desirable  to  make  an  approximate  evaluation  of  the 
parameters  before  using  this  algorithm  in  order  to  avoid  using 
meaningless  values  of  the  parameters. 
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In  their  study  of  the  effect  of  capillary  pressure  and 
gravity  on  two-phase  fluid  flow  in  a  porous  medium,  Fayers  and 
Sheldon  (4)  derived  the  immiscible  displacement  equation  in  its 
Lagrangian  form, 


aX  _  dG 


aT 


dS 


as 


C(S)\ 
3X  I 
3S  / 


(A-l) 


Assumming  that  X(S,T),  C(S)  and  G(S)  are  well  behaved 
functions, they  expanded  them  about  S=0,  which  gives: 


X  =  x  +  S01  E  x .  ( t )  S 
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The  velocity  of  the  foot  of  the  front  is  thus. 


11m  (It)  =  11m  (c-  o+l) 


C,S 
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S->0 
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and,  in  order  for  it  to  be  finite,  the  condition  C=o  is  required. 
Eq.(A-l)  was  then  rewritten  as, 

3X  (3X,2  =  dG  (3X}2  .  c  !fX  +  |C|X  {A-6) 

aT  (asJ  as  3S2  3S  3S 


Combining  Eq.(A-2),  Eq.(A-3),  Eq.(A-4)  and  Eq.(A-6)  yields: 
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With  the  condition  a=c,  Eq.(A-7)  may  be  rewritten  as. 
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Hence,  dividing  both  sides  by  S 


,  and  rearranging, 
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By  comparing  powers  of  S  on  both  sides  of  Eq.(A-9),  and  supposing 
that  a  and  e  are  non-integral,  it  is  easy  to  see  that  either 


e  (e+i-i)  g,  s1'1  =  o 
i-l  1 


or 


a 


The  first  condition  may  be  rejected  because  it  is  equivalent  to 

dX  dG  =  n 

dt  dS  u 

and,  therefore,  the  condition  a  =  c  =  $-1  must  hold.  Fayers  and 
Sheldon  (92),  then,  showed  that,  as  S  goes  to  zero, 

dPc  _  cC-3  dPc  _  o 

dS  ~  *  tnat  1S’  dS  =  S  * 

This  condition,  however,  does  not  originate  from  the  physics  of  the 
immiscible  displacement,  but  rather  from  the  chosen  mathematical 
tool,  the  expansion  in  series  of  X,  C  and  G.  Other  mathematical 
tools,  such  as  finite  differences,  may  be  used  to  solve  the 
immiscible  displacement  equation;  therefore,  it  is  not  necessary  to 
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resort  to  expansions  in  series. 

By  definition  (104) 

,,,,  .  Krw'S>  Krnv/S>  d\ 

{  )  Krw^  +  WS>  dS 


If  Krw(S)  =  6(S)Sm  and  Krnw(S)  =  n(S)(l-S),  (105), 
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then,  with  ^  , 


lim  C(S)  =  — — r4  Sm‘n 


S->0 


alf +  1 
w 


and  the  condition,  lim  C(S)  -  0,  is  still  met  if  n  <  m. 
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(A- 1 0 ) 


Eq . (5-4)  is,  therefore,  justified. 
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The  drainage  capillary  pressure  is  given  by  Eq.(4-9): 


Pc  -  AApw2(rg  -  j)h 


(B-l) 


Therefore, 
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The  other  solution  has  no  physical  meaning. 
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B.l  DERIVATION  OF  EQ . (4-12) 
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Let  h*  be  located  in  the  g —  slice.  The  saturation  of 

the  first  g-1  slices  is,  then,  still  the  maximum  saturation,  S,  =1. 
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Combining  Eq.(B-3)  and  Eq.(B-5)  yields, 
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B.2  DERIVATION  OF  EQ. (4-1 1 ) 
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The  average  saturation  of  any  slice  j  above  the  g —  slice  is, 
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and  with  Eq. (B-3) , 


(B-8) 
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The  imbibition  capillary  pressure  is  given  by  Eq . (4-23 ) 


Pc  =  [°2(p0-  pw>  +  rw  pw  -  ro  p0] 


(C-l) 


with  D  =  rg  -  h.  Therefore, 
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and, 
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In  the  following  derivations,  PQ  is  the  capillary  pressure  at  the 
outer  end  of  the  core,  at  the  end  of  the  imbibition  experiments. 


C-l.  DERIVATION  OF  EQ. (4-24)  and  EQ. (4-25) 
Eq. (4-21 )  can  be  rewritten  as, 
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(C-4) 


Combining  Eq.(C-4)  and  Eq . (4-10)  yields, 
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and,  wi th  Eq. (C-3) , 
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It  follows  that, 
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C. 2  DERIVATION  OF  EQ. (4-26)  AMD  EQ. (4-27) 
Eq . ( 4-22 )  can  be  rewritten  as 
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Combining  Eq.(C-7)  and  Eq . (4-10)  yields, 
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D.l  DERIVATION  OF  EQ. (5-2) 

Let  PI  be  the  mass  of  water  in  the  slice  at  the  end  of  the 
w 

run. 

Mq  be  the  mass  of  oil  in  the  slice  at  the  end  of  the 
run. 

be  the  mass  of  oil  and  water  contained  in  the  slice 
at  the  end  of  the  run. 

Vp  be  the  pore  volume  of  the  slice. 

Pw  and  pQ  be  respectively  the  densities  of  water  and  oil. 


Then, 
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and, 
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Combining  Eq.(D-l)  and  Eq.(D-2)  yields, 
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With  Eq. (D-l ) ,  one  gets 


(D-4) 


This  equation  represents  the  volume  of  water  contained  in  the  slice 
at  the  end  of  the  run;  therefore, 
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ml 
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If  M‘  is  the  mass  of  water  contained  in  the  100%  water  saturated 
w 

slice,  then. 


(0-6) 


Combining  Eq . (D-5 )  and  Eq . ( D-6 )  yields, 
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m; 

w 


V, 


where  Ap  =  pw 


Eq. (D-7)  is  equivalent  to  Eq . (5-2 ) . 


(D-7) 
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D. 2  DERIVATION  OF  EQ. (5-1 ) 

Let  Pa  be  the  atmospheric  pressure 

P  .  and  P  .  be  the  pressures  at  the  inner  end  in  the 
water  and  oil  phases  respectively. 

Poe  be  the  pressure  in  the  oil  phase  at  the  outer  end 
of  the  core. 

H  ,  h  .  r  ,  r  and  r.  are  defined  in  Figure  D-l . 
o  s  o  e  i 

By  definition, 

Poi  ■  pa  +  “2  [(H0  -  hs)(r0  +  ^1)  po]  (D-8) 


and, 


oe 


H 

P  +  a)2  [  H  (r  +  -4)  p  ] 
a  o  o  2  o 


If  we  assume  that  PQ  ,  the  capillary  pressure  at  the  outer  end 
the  core,  is  zero,  then, 


(0-9) 

of 


wi 


=  P 


oe 


to 


CMr„  -  -4-)  pJ 


w- 


and,  with  Eq. (D-9) , 


Pwi  =  Pa  +  “2  [Ho(ro  +  T0  po  -  hs(re  '  pw] 

By  definition,  the  capillary  pressure  at  the  inner  end  of  the 
core  is: 


(D-l 0) 


(D-l 1 ) 


“2  ths^e  - 


+  Ho 


But, 


r-N 


* 
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FIGURE  D-l 


DRAINAGE  CAPILLARY  PRESSURE  CELL 
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therefore, 


(D-12) 


Eq . (D-l 2)  is  equivalent  to  Eq.(5-1) 

D. 3  DERIVATION  OF  EQ. (4-23) 

The  different  terms  used  in  the  following  equations  are 
defined  in  Figure  D-2. 

Let  Pnv/and  Pw  be  the  pressures  in  the  oil  and  water  phases 
respectively,  at  the  point  where  the  capillary  pressure,  P  ,  is  to 
be  calculated.  By  definition, 


(D-l 3) 


(D-l 4) 


By  design,  the  water  and  the  oil  reservoirs  are  opened  to  the 
atmosphere. 
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Free  oil  surface 


^/\  To  the  rotary  axis 
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Free  water  surface 
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FIGURE  D-2  IMBIBITION  CAPILLARY  PRESSURE  CELL 
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The  capillary  pressure  is,  then, 


or 


Wo  +  T>  po  '  hw(rw  +  pw] 


=  (O' 


[(°  -  ro)(ro  +  p0  +  (D  -  rw)(rw  +  4^)  pw] 


D 2  -  r  2  D2  -  r  2 

.  -  ( 
o 


“2  [(  r-2-)  p0  -  (—2-^)  pJ 


w 


=  [d2(p  -  p  )  +  r  2  p  -  r  2  p  ] 

2  L  ^  o  Mw'  w  Hw  o  MoJ 


T  re  -  h)2^0  ‘  PJ  +  rw2  pw  ‘  ro2  po] 


(D-l 5) 


Eq . ( D-l 5)  is  equivalent  to  Eq . (4-23 ) . 


APPENDIX  E 


PARAMETERS  OF  THE  EXPERIMENTS 
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TABLE 

E-l : 

RUN  #1 

Sample 

Sandstone 

Number  of  si  ices  per 

core: 

6 

Speed  of  rotation 

2310  rpm 

for  50  hours 

Length  of  si  ices ,  cm 

#1: 

1.208 

#7:  1.204 

#2: 

1.206 

#8:  1.200 

#3: 

1.213 

#9:  1.207 

#4: 

1.202 

#10:  1.211 

#5: 

1.198 

#11:  1.213 

#6: 

1.206 

#12:  1.220 

TABLE 

E-2: 

RUN  #2 

Sampl e 

P-12-C 

Number  of  si  ices  per 

core: 

3 

Speed  of  rotation 

1035  rpm 

for  15  hours 

1500  rpm 

for  20  hours 

2520  rpm 

for  7  hours 

Length  of  slices,  cm 

#1: 

2.403 

#4:  2.401 

#2: 

2.402 

#5:  2.402 

#3: 

2.405 

#6:  2.407 

2.407 
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TABLE 

E-3: 

RUN  #3 

Sample 

P-12-C 

Number  of  slices  per 

core: 

4 

Speed  of  rotation 

1000  rpm 

for  45 

hours 

1500  rpm 

for  42 

hours 

2500  rpm 

for  46 

hours 

Length  of  si  ices,  cm 

#1: 

2.402 

#5: 

2.403 

#2: 

2.403 

#6: 

2.401 

#3: 

1.900 

#7: 

1.902 

#4: 

0.500 

#8: 

0.498 

TABLE  E-4: 

RUN  #4 

Sample 

P-12-C 

Number  of  slices  per 

core: 

7 

Speed  of  rotation 

2820  rpm 

for  39 

hours 

Length  of  slices,  cm 

#1: 

1.204 

#8: 

1.201 

#2: 

1.209 

#9: 

1.201 

#3: 

1.206 

#10: 

1.204 

#4: 

1.207 

#11: 

1.205 

#5: 

1.200 

#12: 

1.204 

#6: 

0.700 

#13: 

0.701 

#7: 

0.500 

#14: 

0.498 

X 
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TABLE 

E-5: 

RUN  #5 

Sample 

P-12-C 

Number  of  si  ices  per 

core: 

6 

Speed  of  rotation 

2480  rpm 

for  44 

hours 

Length  of  slices,  cm 

#1: 

1.204 

#7: 

1.204 

#2: 

1.209 

#8: 

1.201 

#3: 

1.206 

#9: 

1.201 

#4: 

1.207 

#10: 

1.204 

#5: 

1.200 

#11: 

1.205 

#6: 

1.205 

#12: 

1.204 

TABLE 

E-6: 

RUN  #6 

Sampl e 

P-12-C 

Number  of  slices  per 

core: 

6 

Speed  of  rotation 

2385  rpm 

for  50 

hours 

Length  of  si  ices ,  cm 

#1: 

1.204 

#7: 

1.204 

#2: 

1.209 

#8: 

1.201 

#3: 

1.206 

#9: 

1.201 

#4: 

1.207 

#10: 

1.204 

#5: 

1.200 

#11: 

1.205 

#6: 

1.205 

#12: 

1.204 
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TABLE 

E-7: 

RUN  #7 

Sample 

P-6-C 

Number  of  slices  per 

core: 

3 

Speed  of  rotation 

1360  rpm 

for  20 

hours 

2120  rpm 

for  47 

hours 

2920  rpm 

for  22 

hours 

Length  of  slices,  cm 

#1: 

2.397 

#4: 

2.402 

#2: 

2.385 

#5: 

2.400 

#3: 

2.402 

#6: 

2.403 

TABLE  E-8: 

RUN  #8 

Sampl e 

P-6-C 

Number  of  slices  per 

core: 

6 

Speed  of  rotation 

2910  rpm 

for  216 

hours 

Length  of  slices,  cm 

#1: 

1.195 

#7: 

1.193 

#2: 

1.194 

#8: 

1.191 

#3: 

1.198 

#9: 

1.158 

#4: 

1.199 

#10: 

1.194 

#5: 

1.193 

#11: 

1.193 

#6: 

1.197 

#12: 

1.161 

X 
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EXPERIMENTAL  AND  CALCULATED  DATA 
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TABLE  F-l :  RUN  #1 


h,  cm 

Pc,  kPa 

Experimental 

K 

Calculated 

1.206 

25.55 

0.6165 

0.6050 

1.220 

25.84 

0.5900 

0.6016 

2.405 

49.19 

0.4222 

0.4123 

2.433 

49.73 

0.3982 

0.4094 

3.606 

69.89 

0.3307 

0.3232 

3.644 

71.80 

0.3131 

0.3214 

4.819 

91.50 

0.2837 

0.2747 

4.852 

92.02 

0.2690 

0.2737 

6.025 

109.98 

0.2519 

0.2450 

6.052 

110.36 

0.2391 

0.2445 

7.233 

126.70 

0.2324 

0.2250 

7.256 

127.00 

0.2136 

0.2247 
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TABLE  F-2:  RUN  #2 


Experimental 

Calculated 

h,  cm 

u,  rpm 

pc  kPa 

w 

2.401 

1035 

9.86 

0.9823 

0.9997 

2.403 

1035 

9.87 

0.9828 

0.9996 

4.803 

1035 

18.32 

0.9198 

0.9194 

4.805 

1035 

18.33 

0.9319 

0.9193 

2.401 

1500 

20.71 

0.9036 

0.8944 

2.403 

1500 

20.73 

0.9003 

0.8943 

7.210 

i  035 

25.37 

0.8027 

0.8218 

7.210 

1035 

25.37 

0.8141 

0.8218 

4.803 

1500 

38.48 

0.6871 

0.6918 

4.805 

1500 

38.49 

0.6849 

0.6916 

7.210 

1500 

53.30 

0.5611 

0.5647 

7.210 

1500 

53.30 

0.5717 

0.5647 

2.401 

2520 

58.46 

0.5698 

0.5594 

2.403 

2520 

58.51 

0.5688 

0.5591 

4.803 

2520 

108.60 

0.3769 

0.3763 

4.805 

2520 

108.63 

0.3662 

0.3762 

7.210 

2520 

150.43 

0.2951 

0.3066 

7.210 

2520 

150.44 

0.3198 

0.3066 

1  %  | 

* 

132 


TABLE  F-3:  RUN  #3 


h,  cm 

to,  rpm 

Pc,  kPa 

Experimental 

Cal  culated 

K 

2.400 

1000 

9.20 

0.9926 

0.9965 

2.400 

1000 

9.20 

0.9866 

0.9965 

4.801 

1000 

17.09 

0.9426 

0.9383 

4.803 

1000 

17.10 

0.9491 

0.9383 

2.400 

1500 

20.71 

0.8834 

0.8946 

2.400 

1500 

20.71 

0.8907 

0.8947 

7.204 

1000 

23.67 

0.8255 

0.8393 

4.801 

1500 

38.46 

0.6834 

0.6769 

4.803 

1500 

38.47 

0.6807 

0.6767 

7.204 

1500 

53.26 

0.5604 

0.5463 

2.400 

2500 

57.52 

0.5528 

0.5471 

2.400 

2500 

57.52 

0.5334 

0.5472 

4.801 

2500 

106.84 

0.3571 

0.3694 

4.803 

2500 

106.87 

0.3697 

0.3694 

7.204 

2500 

147.95 

0.2978 

0.3044 

7.205 

2500 

147.97 

0.3174 

0.3044 

0.498 

1000 

2.02 

0.9816 

0.500 

1000 

2.03 

0.9825 

0.498 

1500 

4.54 

0.9461 

0.500 

1500 

4.56 

0.9568 

0.498 

2500 

12.62 

0.8493 

0.500 

2500 

12.66 

0.8736 

< 

'"'X 
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TABLE  F-4: 

RUN  #4 

h,  cm 

P  ,  kPa 

Experimental 

r3 

Cal culated 

r 

c5 

w 

w 

0.498 

16.05 

0.7376 

0.7341 

0.500 

16.11 

0.7319 

0.7335 

1.200 

37.89 

0.5580 

0.5711 

1.200 

37.91 

0.5769 

0.5713 

2.400 

73.18 

0.4040 

0.3945 

2.404 

73.30 

0.3973 

0.3941 

3.606 

106.02 

0.2954 

0.2929 

3.609 

106.11 

0.2834 

0.2928 

4.813 

136.22 

0.2328 

0.2311 

4.814 

136.24 

0.2249 

0.2311 

6.015 

163.67 

0.1882 

0.1912 

6.021 

163.82 

0.1967 

0.1910 

7.216 

188.49 

0.1598 

0.1637 

7.225 

188.68 

0.1688 

0.1635 

V"  X 

TABLE  F-5:  RUN  #5 


h,  cm 

Pc,  kPa 

Experimental 

s7 

Calculated 

K 

1.204 

29.41 

0.7292 

0.7418 

1.205 

29.43 

0.7505 

0.7414 

2.404 

56.70 

0.5404 

0.5347 

2.410 

56.82 

0.5401 

0.5339 

3.611 

82.09 

0.3953 

0.4041 

3.614 

82.15 

0.4040 

0.4039 

4.815 

105.40 

0.3225 

0.3219 

4.817 

105.44 

0.3161 

0.3218 

6.016 

126.61 

0.2736 

0.2687 

6.026 

126.77 

0.2661 

0.2683 

7.220 

145.84 

0.2365 

0.2325 

7.230 

145.99 

0.2309 

0.2322 

< 


rv. 

TABLE  F-6:  RUN  #6 


h,  cm 

Pc,  kPa 

Experimental 

Calculated 

1.204 

27.21 

0.7704 

0.7662 

1.205 

27.21 

0.7582 

0.7657 

2.404 

52.44 

0.5732 

0.5695 

2.410 

52.55 

0.5747 

0.5686 

3.611 

75.93 

0.4360 

0.4363 

3.614 

75.98 

0.4309 

0.4361 

4.815 

97.47 

0.3426 

0.3478 

4.817 

97.52 

0.3482 

0.3477 

6.016 

117.10 

0.2859 

0.2885 

6.026 

117.25 

0.2916 

0.2881 

7.220 

134.88 

0.2470 

0.2477 

7.230 

135.02 

0.2508 

0.2474 

» 
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TABLE  F-7 :  RUN  #7 


h,  cm 

to,  rpm 

PC*  kPa 

Experimental 

r 

w 

Calculated 

!r 

w 

2.397 

1360 

17.00 

0.9906 

0.9755 

2.402 

1360 

17.04 

0.9884 

0.9760 

4.782 

1360 

31.51 

0.9723 

0.9677 

4.802 

1360 

31.62 

0.9768 

0.9670 

2.397 

2120 

41.31 

0.8910 

0.9078 

2.402 

2120 

41.40 

0.8856 

0.9073 

7.184 

1360 

43.69 

0.9115 

0.8780 

7.205 

1360 

43.79 

0.9186 

0.8772 

4.782 

2120 

76.57 

0.6732 

0.6773 

4.802 

2120 

76.84 

0.6802 

0.6758 

2.397 

2920 

78.38 

0.6778 

0.6800 

2.402 

2920 

78.53 

0.6602 

0.6792 

7.184 

2120 

106.18 

0.5292 

0.5358 

7.205 

2120 

106.41 

0.5326 

0.5349 

4.782 

2920 

145.26 

0.4541 

0.4483 

4.802 

2920 

145.78 

0.4319 

0.4472 

7.184 

2920 

201.43 

0.3726 

0.3570 

7.205 

2920 

201.87 

0.3554 

0.3565 
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TABLE  F-8: 

RUN  #8 

h,  cm 

P  ,  kPa 

Experimental 

r 

Cal cul ated 

ir 

c* 

w 

w 

1.161 

39.08 

0.8572 

0.8651 

1.197 

40.25 

0.8686 

0.8567 

2.354 

76.53 

0.6368 

0.6371 

2.389 

77.62 

0.6207 

0.6321 

3.548 

111.26 

0.5135 

0.5027 

3.589 

112.42 

0.4918 

0.4992 

4.706 

142.31 

0.4410 

0.4227 

4.787 

144.40 

0.4100 

0.4183 

5.897 

171.54 

0.3821 

0.3687 

5.981 

173.51 

0.3548 

0.3657 

7.090 

198.08 

0.335* 

0.3314 

7.176 

199.87 

0.3172 

0.3291 

■ 


